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Abstract 

The status of experimental tests of general relativity and of theoretical frameworks for 
analysing them is reviewed. Einstein's equivalence principle (EEP) is well supported by ex- 
periments such as the Eotvos experiment, tests of special relativity, and the gravitational 
redshift experiment. Ongoing tests of EEP and of the inverse square law are searching for 
new interactions arising from unification or quantum gravity. Tests of general relativity at the 
post-Newtonian level have reached high precision, including the light deflection, the Shapiro 
time delay, the perihelion advance of Mercury, and the Nordtvedt effect in lunar motion. 
Gravitational-wave damping has been detected in an amount that agrees with general relativ- 
ity to better than half a percent using the Hulse- Taylor binary pulsar, and other binary pulsar 
systems have yielded other tests, especially of strong-field effects. When direct observation of 
gravitational radiation from astrophysical sources begins, new tests of general relativity will 
be possible. 
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1 Introduction 



At the time of the birth of general relativity (GR), experimental confirmation was almost a side 
issue. Einstein did calculate observable effects of general relativity, such as the perihelion advance 
of Mercury, which he knew to be an unsolved problem, and the deflection of light, which was 
subsequently verified. But compared to the inner consistency and elegance of the theory, he 
regarded such empirical questions as almost peripheral. Today, experimental gravitation is a 
major component of the field, characterized by continuing efforts to test the theory's predictions, 
to search for gravitational imprints of high-energy particle interactions, and to detect gravitational 
waves from astronomical sources. 

The modern history of experimental relativity can be divided roughly into four periods: Genesis, 
Hibernation, a Golden Era, and the Quest for Strong Gravity. The Genesis (1887-1919) comprises 
the period of the two great experiments which were the foundation of relativistic physics - the 
Michelson-Morley experiment and the Eotvos experiment - and the two immediate confirmations 
of GR - the deflection of light and the perihelion advance of Mercury. Following this was a period 
of Hibernation (1920-1960) during which theoretical work temporarily outstripped technology and 
experimental possibilities, and, as a consequence, the field stagnated and was relegated to the 
backwaters of physics and astronomy. 

But beginning around 1960, astronomical discoveries (quasars, pulsars, cosmic background 
radiation) and new experiments pushed GR to the forefront. Experimental gravitation experienced 
a Golden Era (1960-1980) during which a systematic, world-wide effort took place to understand 
the observable predictions of GR, to compare and contrast them with the predictions of alternative 
theories of gravity, and to perform new experiments to test them. The period began with an 
experiment to confirm the gravitational frequency shift of light (1960) and ended with the reported 
decrease in the orbital period of the Hulse- Taylor binary pulsar at a rate consistent with the GR 
prediction of gravity wave energy loss (1979). The results all supported GR, and most alternative 
theories of gravity fell by the wayside (for a popular review, see |282) ~). 

Since 1980, the field has entered what might be termed a Quest for Strong Gravity. Many 
of the remaining interesting weak-field predictions of the theory are extremely small and difficult 
to check, in some cases requiring further technological development to bring them into detectable 
range. The sense of a systematic assault on the weak-field predictions of GR has been supplanted 
to some extent by an opportunistic approach in which novel and unexpected (and sometimes 
inexpensive) tests of gravity have arisen from new theoretical ideas or experimental techniques, 
often from unlikely sources. Examples include the use of laser-cooled atom and ion traps to perform 
ultra-precise tests of special relativity; the proposal of a "fifth" force, which led to a host of new 
tests of the weak equivalence principle; and recent ideas of large extra dimensions, which have 
motived new tests of the inverse square law of gravity at sub-millimeter scales. 

Instead, much of the focus has shifted to experiments which can probe the effects of strong 
gravitational fields. The principal figure of merit that distinguishes strong from weak gravity is 
the quantity e ~ GM/Rc 2 , where G is the Newtonian gravitational constant, M is the characteristic 
mass scale of the phenomenon, R is the characteristic distance scale, and c is the speed of light. 
Near the event horizon of a non-rotating black hole, or for the expanding observable universe, 
e ~ 0.5; for neutron stars, e ~ 0.2. These are the regimes of strong gravity. For the solar system 
e < 10 -5 ; this is the regime of weak gravity. At one extreme are the strong gravitational fields 
associated with Planck-scale physics. Will unification of the forces, or quantization of gravity at 
this scale leave observable effects accessible by experiment? Dramatically improved tests of the 
equivalence principle, of the inverse square law, or of local Lorentz invariancc are being mounted, to 
search for or bound the imprinted effects of Planck-scale phenomena. At the other extreme are the 
strong fields associated with compact objects such as black holes or neutron stars. Astrophysical 
observations and gravitational wave detectors are being planned to explore and test GR in the 
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strong-field, highly-dynamical regime associated with the formation and dynamics of these objects. 

In this Living Review, we shall survey the theoretical frameworks for studying experimental 
gravitation, summarize the current status of experiments, and attempt to chart the future of the 
subject. We shall not provide complete references to early work done in this field but instead 
will refer the reader to the appropriate review articles and monographs, specifically to Theory 
and Experiment in Gravitational Physics |281) . hereafter referred to as TEGP. Additional recent 
reviews in this subject are EHlEHEl E3 El EH • References to TEGP will be by chapter or 
section, e.g. "TEGP 8.9 [2HU" • 
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2 Tests of the Foundations of Gravitation Theory 



2.1 The Einstein equivalence principle 

The principle of equivalence has historically played an important role in the development of gravi- 
tation theory. Newton regarded this principle as such a cornerstone of mechanics that he devoted 
the opening paragraph of the Principia to it. In 1907, Einstein used the principle as a basic ele- 
ment in his development of general relativity. We now regard the principle of equivalence as the 
foundation, not of Newtonian gravity or of GR, but of the broader idea that spacetime is curved. 
Much of this viewpoint can be traced back to Robert Dicke, who contributed crucial ideas about 
the foundations of gravitation theory between 1960 and 1965. These ideas were summarized in 
his influential Lcs Houchcs lectures of 1964 |96| . and resulted in what has come to be called the 
Einstein equivalence principle (EEP). 

One elementary equivalence principle is the kind Newton had in mind when he stated that 
the property of a body called "mass" is proportional to the "weight" , and is known as the weak 
equivalence principle (WEP). An alternative statement of WEP is that the trajectory of a freely 
falling "test" body (one not acted upon by such forces as electromagnetism and too small to be 
affected by tidal gravitational forces) is independent of its internal structure and composition. In 
the simplest case of dropping two different bodies in a gravitational field, WEP states that the 
bodies fall with the same acceleration (this is often termed the Universality of Free Fall, or UFF). 

The Einstein equivalence principle (EEP) is a more powerful and far-reaching concept; it states 
that: 

1. WEP is valid. 

2. The outcome of any local non-gravitational experiment is independent of the velocity of the 
freely-falling reference frame in which it is performed. 

3. The outcome of any local non-gravitational experiment is independent of where and when in 
the universe it is performed. 

The second piece of EEP is called local Lorentz invariance (LLI), and the third piece is called 
local position invariance (LPI). 

For example, a measurement of the electric force between two charged bodies is a local non- 
gravitational experiment; a measurement of the gravitational force between two bodies (Cavendish 
experiment) is not. 

The Einstein equivalence principle is the heart and soul of gravitational theory, for it is pos- 
sible to argue convincingly that if EEP is valid, then gravitation must be a "curved spacetime" 
phenomenon, in other words, the effects of gravity must be equivalent to the effects of living in a 
curved spacetime. As a consequence of this argument, the only theories of gravity that can fully 
embody EEP are those that satisfy the postulates of "metric theories of gravity" , which are: 

1. Spacetime is endowed with a symmetric metric. 

2. The trajectories of freely falling test bodies are geodesies of that metric. 

3. In local freely falling reference frames, the non-gravitational laws of physics are those written 
in the language of special relativity. 

The argument that leads to this conclusion simply notes that, if EEP is valid, then in local 
freely falling frames, the laws governing experiments must be independent of the velocity of the 
frame (local Lorentz invariance), with constant values for the various atomic constants (in order 
to be independent of location). The only laws we know of that fulfill this are those that are 
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compatible with special relativity, such as Maxwell's equations of electromagnetism. Furthermore, 
in local freely falling frames, test bodies appear to be unaccelerated, in other words they move 
on straight lines; but such "locally straight" lines simply correspond to "geodesies" in a curved 
spacetime (TEGP 2.3 [M]). 

General relativity is a metric theory of gravity, but then so are many others, including the 
Brans-Dicke theory and its generalizations. Theories in which varying non-gravitational constants 
are associated with dynamical fields that couple to matter directly are not metric theories. Neither, 
in this narrow sense, is superstring theory (see Sec. 12.31) . which, while based fundamentally on a 
spacetime metric, introduces additional fields (dilatons, moduli) that can couple to material stress- 
energy in a way that can lead to violations, say, of WEP. It is important to point out, however, that 
there is some ambiguity in whether one treats such fields as EEP-violating gravitational fields, or 
simply as additional matter fields, like those that carry electromagnetism or the weak interactions. 
Still, the notion of curved spacetime is a very general and fundamental one, and therefore it is 
important to test the various aspects of the Einstein equivalence principle thoroughly. We first 
survey the experimental tests, and describe some of the theoretical formalisms that have been 
developed to interpret them. For other reviews of EEP and its experimental and theoretical 
significance, see |127lllfi3| . 

2.1.1 Tests of the weak equivalence principle 

A direct test of WEP is the comparison of the acceleration of two laboratory-sized bodies of different 
composition in an external gravitational field. If the principle were violated, then the accelerations 
of different bodies would differ. The simplest way to quantify such possible violations of WEP in 
a form suitable for comparison with experiment is to suppose that for a body with inertial mass 
mi, the passive gravitational mass mp is no longer equal to mi, so that in a gravitational field 
g, the acceleration is given by mia = mpg. Now the inertial mass of a typical laboratory body 
is made up of several types of mass-energy: rest energy, electromagnetic energy, weak-interaction 
energy, and so on. If one of these forms of energy contributes to mp differently than it does to mi, 
a violation of WEP would result. One could then write 



where E A is the internal energy of the body generated by interaction A, and r\ A is a dimensionless 
parameter that measures the strength of the violation of WEP induced by that interaction, and c 
is the speed of light. A measurement or limit on the fractional difference in acceleration between 
two bodies then yields a quantity called the "Eotvos ratio" given by 



where we drop the subscript "I" from the inertial masses. Thus, experimental limits on r\ place 
limits on the WEP-violation parameters t) A . 

Many high-precision Eotvos-type experiments have been performed, from the pendulum exper- 
iments of Newton, Bessel and Potter, to the classic torsion-balance measurements of Eotvos |102j . 
Dicke [57], Braginsky and their collaborators. In the modern torsion-balance experiments, 
two objects of different composition are connected by a rod or placed on a tray and suspended 
in a horizontal orientation by a fine wire. If the gravitational acceleration of the bodies differs, 
and this difference has a component perpendular to the suspension wire, there will be a torque 
induced on the wire, related to the angle between the wire and the direction of the gravitational 
acceleration g. If the entire apparatus is rotated about some direction with angular velocity w, the 




(1) 



A 




5 



torque will be modulated with period 2tt/lu. In the experiments of Eotvos and his collaborators, 
the wire and g were not quite parallel because of the centripetal acceleration on the apparatus due 
to the Earth's rotation; the apparatus was rotated about the direction of the wire. In the Dicke 
and Braginsky experiments, g was that of the Sun, and the rotation of the Earth provided the 
modulation of the torque at a period of 24 hr (TEGP 2.4 (a) 281 ]). Be ginning in the late 1980s, 
numerous experiments were carried out primarily to search for a "fifth force" (see Sec. I2.3.TJ1 . but 
their null results also constituted tests of WEP. In the "free-fall Galileo experiment" performed at 
the University of Colorado, the relative free-fall acceleration of two bodies made of uranium and 
copper was measured using a laser interferomctric technique. The "Eot-Wash" experiments car- 
ried out at the University of Washington used a sophisticated torsion balance tray to compare the 
accelerations of various materials toward local topography on Earth, movable laboratory masses, 
the Sun and the galaxy [2501 121) . and have reached levels of 3 x 10 -13 0]. The resulting upper 
limits on r\ are summarized in Figure^(TEGP 14.1 [281 ; for a bibliography of experiments up to 
1991, see [TU5]1. 

A number of projects are in the development or planning stage to push the bounds on r\ even 
lower. The project MICROSCOPE (MICRO-Satellite a Trainee Compensee pour l'Observation du 
Principe d'Equivalence) is designed to test WEP to 10~ 15 . It is being developed by the French 
space agency CNES for a possible launch in March, 2008, for a one-year mission The drag- 
compensated satellite will be in a Sun-synchronous polar orbit at 700 km altitude, with a payload 
consisting of two differential accelerometers, one with elements made of the same material (plat- 
inum), and another with elements made of different materials (platinum and titanium). 

Another, known as Satellite Test of the Equivalence Principle (STEP) [24B1 , is under consider- 
ation as a possible joint effort of NASA and the European Space Agency (ESA), with the goal of 
a 10 -18 test. STEP would improve upon MICROSCOPE by using cryogenic techniques to reduce 
thermal noise, among other effects. At present, STEP (along with a number of variants, called 
MiniSTEP and QuickSTEP) has not been approved by any agency beyond the level of basic design 
studies or supporting research and development. An alternative concept for a space test of WEP 
is Galileo-Galilei pQ, which uses a rapidly rotating differential accelerometer as its basic element. 
Its goal is a bound on rj at the 10~ 13 level on the ground and 10~ 17 in space. 

2.1.2 Tests of local Lorentz invariance 

Although special relativity itself never benefited from the kind of "crucial" experiments, such as 
the perihelion advance of Mercury and the deflection of light, that contributed so much to the 
initial acceptance of GR and to the fame of Einstein, the steady accumulation of experimental 
support, together with the successful merger of special relativity with quantum mechanics, led to 
its being accepted by mainstream physicists by the late 1920s, ultimately to become part of the 
standard toolkit of every working physicist. This accumulation included 

• the classic Michelson-Morley experiment and its descendents |187l 12381 11421 148) . 

• the Ives-Stillwell, Rossi-Hall and other tests of time-dilation [l37l 12301 1104) . 

• tests of the independence of the speed of light of the velocity of the source, using both binary 
X-ray stellar sources and high-energy pions [4^1 17|. 

• tests of the isotropy of the speed of light 521 12281 1160) 

In addition to these direct experiments, there was the Dirac equation of quantum mechanics 
and its prediction of anti-particles and spin; later would come the stunningly successful relativistic 
theory of quantum electrodynamics. 
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Figure 1: Selected tests of the weak equivalence principle, showing bounds on rj, which measures 
fractional difference in acceleration of different materials or bodies. The free-fall and Eot-Wash 
experiments were originally performed to search for a fifth force ( green region, representing many 
experiments) . The blue band shows evolving bounds on rj for gravitating bodies from lunar laser 
ranging (LLR). 
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In 2005, on the 100th anniversary of the introduction of special relativity, one might ask "what 
is there to test?" . Special relativity has been so thoroughly integrated into the fabric of modern 
physics that its validity is rarely challenged, except by cranks and crackpots. It is ironic then, that 
during the past several years, a vigorous theoretical and experimental effort has been launched, 
on an international scale, to find violations of special relativity. The motivation for this effort is 
not a desire to repudiate Einstein, but to look for evidence of new physics "beyond" Einstein, such 
as apparent violations of Lorentz invariance that might result from certain models of quantum 
gravity. Quantum gravity asserts that there is a fundamental length scale given by the Planck 
length, L p = (hG/c 3 ) 1 ^ 2 — 1.6 x 10~ 33 cm, but since length is not an invariant quantity (Lorentz- 
FitzGerald contraction), then there could be a violation of Lorentz invariance at some level in 
quantum gravity. In brane world scenarios, while physics may be locally Lorentz invariant in 
the higher dimensional world, the confinement of the interactions of normal physics to our four- 
dimensional "brane" could induce apparent Lorentz violating effects. And in models such as string 
theory, the presence of additional scalar, vector and tensor long-range fields that couple to matter 
of the standard model could induce effective violations of Lorentz symmetry. These and other ideas 
have motivated a serious reconsideration of how to test Lorentz invariance with better precision 
and in new ways. 

A simple and useful way of interpreting some of these modern experiments, called the c 2 - 
formalism, is to suppose that the electromagnetic interactions suffer a slight violation of Lorentz 
invariance, through a change in the speed of electromagnetic radiation c relative to the limiting 
speed of material test particles (co, made to take the value unity via a choice of units), in other 
words, c ^ 1 (see Sec. 12.2.31 . Such a violation necessarily selects a preferred universal rest frame, 
presumably that of the cosmic background radiation, through which we are moving at about 
370 km/s |168) . Such a Lorentz-non-invariant electromagnetic interaction would cause shifts in 
the energy levels of atoms and nuclei that depend on the orientation of the quantization axis of 
the state relative to our universal velocity vector, and on the quantum numbers of the state. The 
presence or absence of such energy shifts can be examined by measuring the energy of one such 
state relative to another state that is either unaffected or is affected differently by the supposed 
violation. One way is to look for a shifting of the energy levels of states that are ordinarily equally 
spaced, such as the Zeeman-split 2 J + 1 ground states of a nucleus of total spin J in a magnetic 
field; another is to compare the levels of a complex nucleus with the atomic hyperfine levels of 
a hydrogen maser clock. The magnitude of these "clock anisotropics" would be proportional to 
5 = |c 2 1|. 

The earliest clock anisotropy experiments were the Hughes-Drever experiments, performed in 
the period 1959-60 independently by Hughes and collaborators at Yale University, and by Drever 
at Glasgow University, although their original motivation was somewhat different [1321 199] . The 
Hughes-Drever experiments yielded extremely accurate results, quoted as limits on the parameter 
5 = cr 2 — 1 in Figure|3 Dramatic improvements were made in the 1980s using laser-cooled trapped 
atoms and ions |216l 11641 153] , This technique made it possible to reduce the broading of resonance 
lines caused by collisions, leading to improved bounds on 8 shown in Figure^ (experiments labelled 
NIST, U. Washington and Harvard, respectively). 

Also included for comparison is the corresponding limit obtained from Michelson-Morley type 
experiments (for a review, see [128 ]). In those experiments, when viewed from the preferred frame, 
the speed of light down the two arms of the moving interferometer is c, while it can be shown using 
the electrodynamics of the c 2 formalism, that the compensating Lorentz-FitzGerald contraction of 
the parallel arm is governed by the speed Co = 1. Thus the Michelson-Morley experiment and its 
descendants also measure the coefficient c~ 2 — 1. One of these is the Brillet-Hall experiment |48) . 
which used a Fabry-Perot laser interferometer. In a recent series of experiments, the frequencies 
of electromagnetic cavity oscillators in various orientations were compared with each other or with 
atomic clocks as a function of the orientation of the laboratory |297l 11691 11911 1141 1249) . These 
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placed bounds on c~ 2 — 1 at the level of better than a part in 10 9 . Haugan and Lammerzahl |126| 
have considered the bounds that Michelson-Morley type experiments could place on a modified 
electrodynamics involving a "vector-valued" effective photon mass. 

TESTS OF 
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Figure 2: Selected tests of local Lorentz invariance showing the bounds on the parameter 5, which 
measures the degree of violation of Lorentz invariance in electromagnetism. The Michelson-Morley, 
Joos, Brillet-Hall and cavity experiments test the isotropy of the round-trip speed of light. The 
centrifuge, two-photon absorption ( TP A ) and JPL experiments test the isotropy of light speed 
using one-way propagation. The most precise experiments test isotropy of atomic energy levels. 
The limits assume a speed of Earth of 370 km/s relative to the mean rest frame of the universe. 

The c 2 framework focusses exclusively on classical electrodynamics. It has recently been ex- 
tended to the entire standard model of particle physics by Kostelecky and colleagues [651 1661 H56) . 
The "Standard Model Extension" (SME) has a large number of Lorentz- violating parameters, open- 
ing up many new opportunities for experimental tests (Sec. I2.2.4f> . A variety of clock anisotropy 
experiments have been carried out to bound the electromagnetic parameters of the SME frame- 
work |155| . For example, the cavity experiments described above [2971 11691 1191] placed bounds on 
the coefficients of the tensors k e - and k 0+ (sec Sec. I2.2.4l for definitions) at the levels of 10 -14 and 
10 -10 , respectively. Direct comparisons between atomic clocks based on different nuclear species 
place bounds on SME parameters in the neutron and proton sectors, depending on the nature of 
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the transitions involved. The bounds achieved range from 10~ 27 to 10 -32 GeV. 

Astrophysical observations have also been used to bound Lorentz violations. For example, if 
photons satisfy the Lorentz violating dispersion relation 

E 2 = pV + E Pl fW \p\c + /(Vc 2 + W 3 C 3 + . . . , (3) 

where E Pl = {hc 5 /G) 1/2 is the Planck energy, then the speed of light v 7 = dE/dp would be given, 
to linear order in the by 

w 7 ^(n-l)& ) E n - 2 

- 1 « 1 + > . 4 

n>l * L 

Such a Lorentz-violating dispersion relation could be a relic of quantum gravity, for instance. 
By bounding the difference in arrival time of high-energy photons from a burst source at large 
distances, one could bound contributions to the dispersion for n > 2. One limit, |/' 3 '| < 128 
comes from observations of 1 and 2 TeV gamma rays from the blazar Markarian 421 |22J. Another 
limit comes from birefringence in photon propagation: in many Lorentz violating models, different 
photon polarizations may propagate with different speeds, causing the plane of polarization of 
a wave to rotate. If the frequency dependence of this rotation has a dispersion relation similar 
to Eq. then by studying "polarization diffusion" of light from a polarized source in a given 
bandwidth, one can effectively place a bound |/^| < 10~ 4 |120| . Other testable effects of Lorentz 
invariance violation include threshold effects in particle reactions, gravitational Ccrenkov radiation, 
and neutrino oscillations. 

Mattingly |183| gives a thorough and up-to-date review of both the theoretical frameworks and 
the experimental results for tests of LLI. 



2.1.3 Tests of local position invariance 

The principle of local position invariance, the third part of EEP, can be tested by the gravitational 
redshift experiment, the first experimental test of gravitation proposed by Einstein. Despite the 
fact that Einstein regarded this as a crucial test of GR, we now realize that it does not distinguish 
between GR and any other metric theory of gravity, but is only a test of EEP. A typical gravitational 
redshift experiment measures the frequency or wavelength shift Z = Av/v = — AX/X between two 
identical frequency standards (clocks) placed at rest at different heights in a static gravitational 
field. If the frequency of a given type of atomic clock is the same when measured in a local, 
momentarily comoving freely falling frame (Lorentz frame), independent of the location or velocity 
of that frame, then the comparison of frequencies of two clocks at rest at different locations boils 
down to a comparison of the velocities of two local Lorentz frames, one at rest with respect to one 
clock at the moment of emission of its signal, the other at rest with respect to the other clock at 
the moment of reception of the signal. The frequency shift is then a consequence of the first-order 
Doppler shift between the frames. The structure of the clock plays no role whatsoever. The result 
is a shift 

Z = AU/c 2 , (5) 

where AU is the difference in the Newtonian gravitational potential between the receiver and the 
emitter. If LPI is not valid, then it turns out that the shift can be written 

Z = (1 + a)AU/c 2 , (6) 

where the parameter a may depend upon the nature of the clock whose shift is being measured 
(see TEGP 2.4 (c) 281 for details). 
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The first successful, high-precision redshift measurement was the series of Pound-Rebka-Snider 
experiments of 1960-1965 that measured the frequency shift of gamma-ray photons from 57 Fe as 
they ascended or descended the Jefferson Physical Laboratory tower at Harvard University. The 
high accuracy achieved - one percent - was obtained by making use of the Mossbauer effect to 
produce a narrow resonance line whose shift could be accurately determined. Other experiments 
since 1960 measured the shift of spectral lines in the Sun's gravitational field and the change in 
rate of atomic clocks transported aloft on aircraft, rockets and satellites. Figure |31 summarizes the 
important redshift experiments that have been performed since 1960 (TEGP 2.4 (c) |281j ). 
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Figure 3: Selected tests of local position invariance via gravitational redshift experiments, showing 
bounds on a, which measures degree of deviation of redshift from the formula Av jv = AU/c 2 . In 
null redshift experiments, the bound is on the difference in a between different kinds of clocks. 



After almost 50 years of inconclusive or contradictory measurements, the gravitational redshift 
of solar spectral lines was finally measured reliably. During the early years of GR, the failure 
to measure this effect in solar lines was siezed upon by some as reason to doubt the theory. 
Unfortunately, the measurement is not simple. Solar spectral lines are subject to the "limb effect" , 
a variation of spectral line wavelengths between the center of the solar disk and its edge or "limb" ; 
this effect is actually a Doppler shift caused by complex convective and turbulent motions in the 
photosphere and lower chromosphere, and is expected to be minimized by observing at the solar 
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limb, where the motions are predominantly transverse. The secret is to use strong, symmetrical 
lines, leading to unambiguous wavelength measurements. Successful measurements were finally 
made in 1962 and 1972 (TEGP 2.4 (c) 281 ). In 1991, LoPresto et al. [1731 measured the solar 
shift in agreement with LPI to about 2 percent by observing the oxygen triplet lines both in 
absorption in the limb and in emission just off the limb. 

The most precise standard redshift test to date was the Vessot-Levine rocket experiment that 
took place in June 1976 |264j . A hydrogen-maser clock was flown on a rocket to an altitude of 
about 10,000 km and its frequency compared to a similar clock on the ground. The experiment 
took advantage of the masers' frequency stability by monitoring the frequency shift as a function of 
altitude. A sophisticated data acquisition scheme accurately eliminated all effects of the first-order 
Doppler shift due to the rocket's motion, while tracking data were used to determine the payload's 
location and the velocity (to evaluate the potential difference AU, and the special relativistic time 
dilation). Analysis of the data yielded a limit \a\ < 2 x 10~ 4 . 

A "null" redshift experiment performed in 1978 tested whether the relative rates of two dif- 
ferent clocks depended upon position. Two hydrogen maser clocks and an ensemble of three 
superconducting-cavity stabilized oscillator (SCSO) clocks were compared over a 10-day period. 
During the period of the experiment, the solar potential U/c 2 changed sinusoidally with a 24-hour 
period by 3 x 10 -13 because of the Earth's rotation, and changed linearly at 3 x 10~ 12 per day 
because the Earth is 90 degrees from perihelion in April. However, analysis of the data revealed no 
variations of either type within experimental errors, leading to a limit on the LPI violation param- 
eter |q h — a scso | < 2 x 10~ 2 259 . This bound has been improved using more stable frequency 
standards, such as atomic fountain clocks |121l 12171 I25| . The current bound, from comparing a 
Cesium atomic fountain with a Hydrogen maser for a year, is |a H - a Cs | < 2.1 x 10" 5 [23- 

The varying gravitational redshift of Earth-bound clocks relative to the highly stable millisecond 
pulsar PSR 1937+21, caused by the Earth's motion in the solar gravitational field around the 
Earth-Moon center of mass (amplitude 4000 km), was measured to about 10 percent |252) . Two 
measurements of the redshift using stable oscillator clocks on spacecraft were made at the one 
percent level: one used the Voyager spacecraft in Saturn's gravitational field |159| . while another 
used the Galileo spacecraft in the Sun's field |161| . 

The gravitational redshift could be improved to the 10~ 10 level using an array of laser cooled 
atomic clocks on board a spacecraft which would travel to within four solar radii of the Sun |181) . 

Modern advances in navigation using Earth-orbiting atomic clocks and accurate time-transfer 
must routinely take gravitational redshift and time-dilation effects into account. For example, the 
Global Positioning System (GPS) provides absolute positional accuracies of around 15 m (even 
better in its military mode), and 50 nanoseconds in time transfer accuracy, anywhere on Earth. 
Yet the difference in rate between satellite and ground clocks as a result of relativistic effects is 
a whopping 39 microseconds per day (46 /zs from the gravitational redshift, and — 7 fj,s from time 
dilation). If these effects were not accurately accounted for, GPS would fail to function at its 
stated accuracy. This represents a welcome practical application of GR! (For the role of GR in 
GPS, see [T7J [T5J; for a popular essay, see [387). ) 

Local position invariance also refers to position in time. If LPI is satisfied, the fundamental 
constants of non-gravitational physics should be constants in time. Tabled shows current bounds 
on cosmological variations in selected dimcnsionless constants. For discussion and references to 
early work, see TEGP 2.4 (c) |281| or jlOOj . For a comprehensive recent review both of experiments 
and of theoretical ideas that underly proposals for varying constants, see |262) . 

Experimental bounds on varying constants come in two types: bounds on the present rate of 
variation, and bounds on the difference between today's value and a value in the distant past. 
The main example of the former type is the clock comparison test, in which highly stable atomic 
clocks of different fundamental type are intercompared over periods ranging from months to years 
(variants of the null redshift experiment). If the frequencies of the clocks depend differently on 
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Limit on k/k 








Constant k 


(yr- 1 ) 


Redshift 


Method 


Fine structure 


< 30 x 10- 1H 







Clock comparisons 182 33 112, 210 


constant 


< 0.5 x 10~ 16 


0.15 




Oklo Natural Reactor [TBI ITTTI |2"TT] 


oem = e 2 /hc 


< 3.4 x 10~ 16 


0.45 




187 Re decay in meteorites [206 




(6.4 ± 1.4) x 10~ 16 


0.2- 


3.7 


Spectra in distant quasars |2(>!)| |194| 




< 1.2 x 10" 16 


0.4- 


2.3 


Spectra in distant quasars 12431 1531 


Weak interaction 


< 1 x 10" 11 


0.15 




Oklo Natural Reactor 75 


constant 


< 5 x 10~ 12 


10 9 




Big Bang nucleosynthesis 180 224 


aw = GfWpc/fi, 










e-p mass ratio 


< 3 x 10~ ib 


2.6- 


3.0 


Spectra in distant quasars 136 



Table 1: Bounds on cosmological variation of fundamental constants of non- gravitational physics. 
For an in-depth review, see '2621. 



the electromagnetic fine structure constant Oem, the electron-proton mass ratio m e /m p , or the 
gyromagnetic ratio of the proton g pi for example, then a limit on a drift of the fractional frequency 
difference translates into a limit on a drift of the constant (s). The dependence of the frequencies on 
the constants may be quite complex, depending on the atomic species involved. The most recent 
experiments have exploited the techniques of laser cooling and trapping, and of atom fountains, in 
order to achieve extreme clock stability, and compared the Rubidium-87 hyperfine transition |182) . 
the Mercury-199 ion electric quadrupole transition the atomic Hydrogen IS — 2S transi- 
tion |112) . or an optical transition in Yttcrbium-171 |210| . against the ground-state hyperfine 
transition in Cesium-133. These experiments show that, today, qem/ cvem < 3 x 10 yr . 

The second type of bound involves measuring the relics of or signal from a process that occurred 
in the distant past and comparing the inferred value of the constant with the value measured in 
the laboratory today. One sub-type uses astronomical measurements of spectral lines at large 
redshift, while the other uses fossils of nuclear processes on Earth to infer values of constants early 
in geological history. 

Earlier comparisons of spectral lines of different atoms or transitions in distant galaxies and 
quasars produced bounds «em or g p {m e /m p ) on the order of a part in 10 per Hubble time [298 . 
Dramatic improvements in the precision of astronomical and laboratory spectroscopy, in the ability 
to model the complex astronomical environments where emission and absorption lines are produced, 
and in the ability to reach large redshift have made it possible to improve the bounds significantly. 
In fact, in 1999, Webb et al. [2691 1194] announced that measurements of absorption lines in Mg, 
Al, Si, Cr, Fe, Ni and Zn in quasars in the redshift range 0.5 < Z < 3.5 indicated a smaller 
value of aEM in earlier epochs, namely Ac^em/c^em = (— 0.72 ± 0.18) x 10" 5 , corresponding to 
c^em/ Q!em = (6.4 ± 1.4) x 10 _16 yr _1 (assuming a linear drift with time). Measurements by other 
groups have so far failed to confirm this non-zero effect [213 ESI E20|; a recent analysis of Mg 
absorption systems in quasars at 0.4 < Z < 2.3 gave Aem/ c<em = (—0.6 ± 0.6) x 10~ 16 yr -1 |243| . 

Another important set of bounds arises from studies of the "Oklo" phenomenon, a group of 
natural, sustained 235 U fission reactors that occurred in the Oklo region of Gabon, Africa, around 
1.8 billion years ago. Measurements of ore samples yielded an abnormally low value for the ratio of 
two isotopes of Samarium, 149 Sm/ 147 Sm. Neither of these isotopes is a fission product, but 149 Sm 
can be depleted by a flux of neutrons. Estimates of the neutron fluence (integrated dose) during 
the reactors' "on" phase, combined with the measured abundance anomaly, yield a value for the 
neutron cross-section for 149 Sm 1.8 billion years ago that agrees with the modern value. However, 
the capture cross-section is extremely sensitive to the energy of a low-lying level (E ~ 0.1 eV), so 
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that a variation in the energy of this level of only 20 milli-eV over a billion years would change 
the capture cross-section from its present value by more than the observed amount. This was first 
analysed in 1976 by Shlyakter [233] . Recent reanalyses of the Oklo data [751ITT7IOTT] lead to a 
bound on «em at the level of around 5 x 10 _17 yr —1 . 

In a similar manner, recent re-analyses of decay rates of 187 Re in ancient meteorites (4.5 billion 
years old) gave the bound ctem/ ctem < 3.4 x 10~ 16 yr -1 206 . 

2.2 Theoretical Frameworks for Analyzing EEP 
2.2.1 Schiff's conjecture 

Because the three parts of the Einstein equivalence principle discussed above are so very different in 
their empirical consequences, it is tempting to regard them as independent theoretical principles. 
On the other hand, any complete and self-consistent gravitation theory must possess sufficient 
mathematical machinery to make predictions for the outcomes of experiments that test each prin- 
ciple, and because there are limits to the number of ways that gravitation can be meshed with the 
special relativistic laws of physics, one might not be surprised if there were theoretical connections 
between the three sub-principles. For instance, the same mathematical formalism that produces 
equations describing the free fall of a hydrogen atom must also produce equations that determine 
the energy levels of hydrogen in a gravitational field, and thereby the ticking rate of a hydrogen 
maser clock. Hence a violation of EEP in the fundamental machinery of a theory that manifests 
itself as a violation of WEP might also be expected to show up as a violation of local position 
invariance. Around 1960, Schiff conjectured that this kind of connection was a necessary feature of 
any self-consistent theory of gravity. More precisely, Schiff's conjecture states that any complete, 
self- consistent theory of gravity that embodies WEP necessarily embodies EEP. In other words, the 
validity of WEP alone guarantees the validity of local Lorentz and position invariance, and thereby 
of EEP. 

If Schiff's conjecture is correct, then Eotvos experiments may be seen as the direct empirical 
foundation for EEP, hence for the interpretation of gravity as a curved-spacetime phenomenon. Of 
course, a rigorous proof of such a conjecture is impossible (indeed, some special counter-examples 
are known [2051 11951 164| ). yet a number of powerful "plausibility" arguments can be formulated. 

The most general and elegant of these arguments is based upon the assumption of energy 
conservation. This assumption allows one to perform very simple cyclic gedanken experiments 
in which the energy at the end of the cycle must equal that at the beginning of the cycle. This 
approach was pioneered by Dicke, Nordtvedt and Haugan (see, e.g. |125| ). A system in a quantum 
state A decays to state B, emitting a quantum of frequency v. The quantum falls a height H in 
an external gravitational field and is shifted to frequency z/, while the system in state B falls with 
acceleration g B . At the bottom, state A is rebuilt out of state B, the quantum of frequency z/, 
and the kinetic energy m B g B H that state B has gained during its fall. The energy left over must 
be exactly enough, niAgAH, to raise state A to its original location. (Here an assumption of local 
Lorentz invariance permits the inertial masses itia and tub to be identified with the total energies 
of the bodies.) If gA and g B depend on that portion of the internal energy of the states that was 
involved in the quantum transition from A to B according to 

g A = .9(1 + aEA/mAC 2 ), g B = g(l + aE B /m B c 2 ), E A - E B = hv (7) 

(violation of WEP), then by conservation of energy, there must be a corresponding violation of 
LPI in the frequency shift of the form (to lowest order in hv/mc 2 ) 

Z = (1/ - v)/v' = (1 + a)gH/c 2 = (1 + a)AU/c 2 . (8) 

Haugan generalized this approach to include violations of LLI 125 , (TEGP 2.5 281 .). 
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Box 1. The THe/j, formalism 



1. Coordinate system and conventions: x = t = time coordinate associated with the 
static nature of the static spherically symmetric (SSS) gravitational field; x = (x,y, z) = 
isotropic quasi-Cartesian spatial coordinates; spatial vector and gradient operations as in 
Cartesian space. 

2. Matter and field variables: 

• rriQa = rest mass of particle a. 

• e a = charge of particle a. 

• x%(t) = world line of particle a. 

• Va — dx£/dt — coordinate velocity of particle a. 

• = electromagnetic vector potential; E = VA — dA/dt, B = VxA 

3. Gravitational potential: f7(x) 

4. Arbitrary functions: T(U), H(U), e(U), fi(U); EEP is satisfied if e = n = (H/T) 1 / 2 for 
all U. 

5. Action: 

/ = -£mo« [(T-HvD^dt + ^ea I A^aZWdt 

a J a J 

+(8tt)- 1 J \eE 2 - ii- x B 2 )d^x. 

6. Non-Metric parameters: 

T Q = -c 2 (d/dU)He(T/H) 1 / 2 ] , 
A = -c 2 (d/dU)H^T/H) 1 / 2 ] , 
T = 1-(TH- 1 €^) , 

where co = (To/ 'Ho) 1 / 2 and subscript "0" refers to a chosen point in space. If EEP is satisfied, 
To = A = T = 0. 
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2.2.2 The THe^i formalism 



The first successful attempt to prove Schiff's conjecture more formally was made by Lightman and 
Lee |167| . They developed a framework called the THe/i formalism that encompasses all metric 
theories of gravity and many non-metric theories (Box^). It restricts attention to the behavior 
of charged particles (electromagnetic interactions only) in an external static spherically symmetric 
(SSS) gravitational field, described by a potential U. It characterizes the motion of the charged 
particles in the external potential by two arbitrary functions T(U) and H(U), and characterizes 
the response of electromagnetic fields to the external potential (gravitationally modified Maxwell 
equations) by two functions e(U) and fJ,(U). The forms of T, H, e and fi vary from theory to 
theory, but every metric theory satisfies 

e = /x = {H/Tf'\ (9) 

for all U. This consequence follows from the action of electrodynamics with a "minimal" or metric 
coupling: 

"lib / ^9^ a 9 v0 F^F aP d A x, (10) 

where the variables are defined in Box^ and where F^ v = A„. M — A^ v . By identifying goo = T and 
gij = HSij in a SSS field, Fjo = -E* and Fij = e^^B^, one obtains Eq. Conversely, every theory 
within this class that satisfies Eq. @ can have its electrodynamic equations cast into "metric" 
form. In a given non-metric theory, the functions T, H, e and fj, will depend in general on the 
full gravitational environment, including the potential of the Earth, Sun and Galaxy, as well as 
on cosmological boundary conditions. Which of these factors has the most influence on a given 
experiment will depend on the nature of the experiment. 

Lightman and Lee then calculated explicitly the rate of fall of a "test" body made up of inter- 
acting charged particles, and found that the rate was independent of the internal electromagnetic 
structure of the body (WEP) if and only if Eq. © was satisfied. In other words, WEP EEP 
and Schiff's conjecture was verified, at least within the restrictions built into the formalism. 

Certain combinations of the functions T, H, e and /i reflect different aspects of EEP. For 
instance, position or [/-dependence of either of the combinations e(T/ ' H) 1 / 2 and n(T/ ' H) 1 / 2 signals 
violations of LPI, the first combination playing the role of the locally measured electric charge or 
fine structure constant. The "non- metric parameters" To and Ao (Box^l are measures of such 
violations of EEP. Similarly, if the parameter To = 1 — (TiJ _1 e/i)o is non-zero anywhere, then 
violations of LLI will occur. This parameter is related to the difference between the speed of light 
c, and the limiting speed of material test particles cq, given by 

c=(eoMor 1/2 , c = (T /Ho) 1 / 2 . (11) 

In many applications, by suitable definition of units, Co can be set equal to unity. If EEP is valid, 
To = Ao = To = everywhere. 

The rate of fall of a composite spherical test body of electromagnetically interacting particles 
then has the form 

a=(mp/m)VU, (12) 



mp/m = 1 + (£f S /Mc 2 ) 



2To 



3 To 



[E^/Mc 2 ) 



2A n 



4 

3 To 



(13) 
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where E^ and are the electrostatic and magnetostatic binding energies of the body, given 
by 




where r a & = |x a , — x&|, n a b = (x a — Xb)/r a b, and the angle brackets denote an expectation value 
of the enclosed operator for the system's internal state. Eotvos experiments place limits on the 
WEP- violating terms in Eq. (|13() . and ultimately place limits on the non- metric parameters |To| < 
2 x 10~ 10 and |Ao| < 3 x 10 -6 . (We set To = because of very tight constraints on it from tests 
of LLI; see Figure [21 where 6 = — T.) These limits are sufficiently tight to rule out a number of 
non-metric theories of gravity thought previously to be viable (TEGP 2.6 (f) |281) L 

The THefi formalism also yields a gravitationally modified Dirac equation that can be used 
to determine the gravitational redshift experienced by a variety of atomic clocks. For the redshift 
parameter a (Eq. ©), the results are (TEGP 2.6 (c) [55T] ): 

{— 3Fo + Ao hydrogen hyperfme transition, H-Maser clock 
— ^(3Fo + Ao) electromagnetic mode in cavity, SCSO clock, (16) 
— 2Po phonon mode in solid, principal transition in hydrogen. 

The redshift is the standard one (a = 0), independently of the nature of the clock if and only 
if Tq = Ao = 0. Thus the Vessot-Levine rocket redshift experiment sets a limit on the parameter 
combination |3Po — Aq| (Figure OJ; the null-redshift experiment comparing hydrogen-maser and 
SCSO clocks sets a limit on |a H - a SC so| = §|r - A Q |. Alvarez and Mann |31 151 ITU1 ITT1 IT2] 
extended the THefi formalism to permit analysis of such effects as the Lamb shift, anomalous 
magnetic moments and non-baryonic effects, and placed interesting bounds on EEP violations. 

2.2.3 The c 2 formalism 

The THefi, formalism can also be applied to tests of local Lorentz invariance, but in this context 
it can be simplified. Since most such tests do not concern themselves with the spatial variation of 
the functions T, H, e, and fi, but rather with observations made in moving frames, we can treat 
them as spatial constants. Then by rescaling the time and space coordinates, the charges and the 
electromagnetic fields, we can put the action in Box^into the form (TEGP 2.6 (a) |281p 

1 = - E m ^ A 1 - <) 1/2dt + E e - I MO<dt + I (E 2 ~ c 2 B 2 )d 4 Xl (17) 

a a J 

where c 2 = iJo/ToeoMo = (1 ~ To)" 1 . This amounts to using units in which the limiting speed Co 
of massive test particles is unity, and the speed of light is c. If c ^ 1, LLI is violated; furthermore, 
the form of the action above must be assumed to be valid only in some preferred universal rest 
frame. The natural candidate for such a frame is the rest frame of the microwave background. 

The electrodynamical equations which follow from Eq. (|17|) yield the behavior of rods and 
clocks, just as in the full THefj, formalism. For example, the length of a rod which moves with 
velocity V relative to the rest frame in a direction parallel to its length will be observed by a 
rest observer to be contracted relative to an identical rod perpendicular to the motion by a factor 
1 — V 2 /2 + 0(V 4 ). Notice that c does not appear in this expression, because only electrostatic 
interactions are involved, and c appears only in the magnetic sector of the action (|17|) . The energy 
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and momentum of an electromagnetically bound body moving with velocity V relative to the rest 
frame are given by 

E = M R + ±M R F 2 + UM^VWi + 0(MV 4 ) , 

221 , 18 , 

P* = M R V l + SM^V 1 + 0(MV 3 ) , 

where Mr = Mo — E^ s , Mo is the sum of the particle rest masses, i?f s is the electrostatic binding 
energy of the system (Eq. (fH| with T^^Hqc^ 1 = 1), and 



where 




Note that (cr 2 — 1) corresponds to the parameter 8 plotted in Figure |21 

The electrodynamics given by Eq. (|17|) can also be quantized, so that we may treat the in- 
teraction of photons with atoms via perturbation theory. The energy of a photon is % times its 
frequency uj, while its momentum is fuv/c. Using this approach, one finds that the difference in 
round trip travel times of light along the two arms of the interferometer in the Michelson-Morley 
experiment is given by Lq(v 2 /c)(c~ 2 — 1). The experimental null result then leads to the bound on 
[cT 2 — 1) shown on Figure[3] Similarly the anisotropy in energy levels is clearly illustrated by the 
tensorial terms in Eqs. (|18|l and l|2U[l : by evaluating E B 13 for each nucleus in the various Hughes- 
Drever-type experiments and comparing with the experimental limits on energy differences, one 
obtains the extremely tight bounds also shown on Figure |21 

The behavior of moving atomic clocks can also be analysed in detail, and bounds on (cr 2 — 1) 
can be placed using results from tests of time dilation and of the propagation of light. In some 
cases, it is advantageous to combine the c 2 framework with a "kinematical" viewpoint that treats a 
general class of boost transformations between moving frames. Such kinematical approaches have 
been discussed by Robertson, Mansouri and Sexl, and Will (see |279p . 

For example, in the "JPL" experiment, in which the phases of two hydrogen masers connected 
by a fiberoptic link were compared as a function of the Earth's orientation, the predicted phase 
difference as a function of direction is, to first order in V, the velocity of the Earth through the 
cosmic background, 

A0/^ » -~(1 - c 2 )(V • n V n ), (21) 

where (f> = 2t:vL 1 v is the maser frequency, L = 21 km is the baseline, and where n and no are 
unit vectors along the direction of propagation of the light at a given time and at the initial time 
of the experiment, respectively. The observed limit on a diurnal variation in the relative phase 
resulted in the bound |c~ 2 — 1| < 3 x 10~ 4 . Tighter bounds were obtained from a "two-photon 
absorption" (TPA) experiment, and a 1960s series of "Mossbauer-rotor" experiments, which tested 
the isotropy of time dilation between a gamma ray emitter on the rim of a rotating disk and an 
absorber placed at the center [279 . 

2.2.4 The Standard Model Extension (SME) 

Kostelecky and collaborators developed a useful and elegant framework for discussing violations 
of Lorentz symmetry in the context of the standard model of particle physics |65l 1661 1156] . Called 
the Standard Model Extension (SME), it takes the standard SU(3) x SU(2) x U(l) field theory of 
particle physics, and modifies the terms in the action by inserting a variety of tensorial quantities 
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in the quark, lepton, Higgs, and gauge boson sectors that could explicitly violate LLI. SME extends 
the earlier classical THe^i and c 2 frameworks, and the x ~ 9 framework of Ni |195| to quantum 
field theory and particle physics. The modified terms split naturally into those that are odd under 
CPT (i.e. that violate CPT) and terms that are even under CPT. The result is a rich and complex 
framework, with many parameters to be analysed and tested by experiment. Such details are 
beyond the scope of this review; for a review of SME and other frameworks, the reader is referred 
to the Living Review by Mattingly |183| . 

Here we confine our attention to the electromagnetic sector, in order to link the SME with the 
c 2 framework discussed above. In the SME, the Lagrangian for a scalar particle 4> with charge e 
interacting with electrodynamics takes the form 

C = [rT + (k^ v ](D^D u <j>-m 2 ^(l) 

~[rTv vP + {kpT^F^F^ , (22) 

where D^cj) = d^<j) + ieA^cf), and where {k^)^ v is a real symmetric trace-free tensor, and (k F )^ a P 
is a tensor with the symmetries of the Riemann tensor, and with vanishing double trace. It has 19 
independent components. There could also be a CPT-odd term in C of the form (kA^c^apA^F 01 ", 
but because of a variety of pre-existing theoretical and experimental constraints, it is generally set 
to zero. 

The tensor (k F )^ au P can be decomposed into "electric", "magnetic" and "odd-parity" compo- 
nents, by defining 

( KDE y k = -2(k F )°^ k , 

(K HB y k = -e jpq e krs (k F ) pqrs , 

{K DB ) kj = -{k HE ) jk =e' pq {k F f kpq . (23) 
In many applications it is useful to use the further decomposition 



K>tr 


= l^ DE y\ 


(Ke+y k 


i, 

= ^{kde + 


KHB) 3k , 


(~Ke-y k 


1, 

= ^(KDE- 


K HB y k - 


(Ho+y k 


l, 

= ^(KDB + 


KHE) 3k , 


(~K -y k 


1, 


K HE ) ik ■ 



(24) 

The first expression is a single number, the next three are symmetric trace-free matrices, and the 
final is an antisymmetric matrix, accounting thereby for the 19 components of the original tensor 

{k F y avf} . 

In the rest frame of the universe, these tensors have some form that is established by the global 
nature of the solutions of the overarching theory being used. In a frame that is moving relative to 
the universe, the tensors will have components that depend on the velocity of the frame, and on 
the orientation of the frame relative to that velocity. 

In the case where the theory is rotationally symmetric in the preferred frame, the tensors (fc</,) A " y 
and (hp) 111 ' can be expressed in the form 

(fc r = Ki(u»u v + ~vr) , 

{k F Y uaP = ^(WV^ufl - if [a rfi v ) , (25) 
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where [] around indices denote antisymmetrization, and where u M is the four- velocity of an observer 
at rest in the preferred frame. With this assumption, all the tensorial quantities in Eq. (|24[) vanish 
in the preferred frame, and, after suitable rescalings of coordinates and fields, the action (|22|l can 
be put into the form of the c 2 framework, with 



1 



1/2 



1 — «tr 

1 + k tr 



1/2 



(26) 



2.3 EEP, particle physics, and the search for new interactions 

Thus far, we have discussed EEP as a principle that strictly divides the world into metric and non- 
metric theories, and have implied that a failure of EEP might invalidate metric theories (and thus 
general relativity) . On the other hand, there is mounting theoretical evidence to suggest that EEP 
is likely to be violated at some level, whether by quantum gravity effects, by effects arising from 
string theory, or by hitherto undetected interactions. Roughly speaking, in addition to the pure 
Einstcinian gravitational interaction, which respects EEP, theories such as string theory predict 
other interactions which do not. In string theory, for example, the existence of such EEP-violating 
fields is assured, but the theory is not yet mature enough to enable a robust calculation of their 
strength relative to gravity, or a determination of whether they are long range, like gravity, or 
short range, like the nuclear and weak interactions, and thus too short range to be detectable. 

In one simple example |95) . one can write the Lagrangian for the low-energy limit of a string- 
inspired theory in the so-called "Einstein frame" , in which the gravitational Lagrangian is purely 
general relativistic: 



C 



U (<p)g' iV 9 al3 F lia F u p 



1> 



(27) 



where g^ u is the non-physical metric, is the Ricci tensor derived from it, ip is a dilaton field, and 
G, U and M are functions of p. The Lagrangian includes that for the electromagnetic field , and 
that for particles, written in terms of Dirac spinors ip. This is not a metric representation because 
of the coupling of ip to matter via M(ip) and U(tp). A conformal transformation g^ v — F(ip)g^ u , 
•tj) = F(ip)~ 3 / 4 ip, puts the Lagrangian in the form ("Jordan" frame) 



A T-i / , d^Fd„F 



-U{p)g^g^F m F K 







ie^^ + ^+qA^-M^F 1 ' 2 



1> 



(28) 



One may choose F(p) — const. /M(ip) 2 so that the particle Lagrangian takes the metric form 
(no explicit coupling to ip) , but the electromagnetic Lagrangian will still couple non-metrically to 
U(<p). The gravitational Lagrangian here takes the form of a scalar-tensor theory (Sec. l3.3^|l . But 
the non-metric electromagnetic term will, in general, produce violations of EEP. For examples of 
specific models, see [2551 188) . Another class of non- metric theories are included in the "varying 
speed of light (VSL)" theories; for a detailed review, see |179) 

On the other hand, whether one views such effects as a violation of EEP or as effects arising from 
additional "matter" fields whose interactions, like those of the electromagnetic field, do not fully 
embody EEP, is to some degree a matter of semantics. Unlike the fields of the standard model of 
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electromagnetic, weak and strong interactions, which couple to properties other than mass-energy 
and are either short range or are strongly screened, the fields inspired by string theory could be 
long range (if they remain massless by virtue of a symmetry, or at best, acquire a very small mass), 
and can couple to mass-energy, and thus can mimic gravitational fields. Still, there appears to be 
no way to make this precise. 

As a result, EEP and related tests are now viewed as ways to discover or place constraints 
on new physical interactions, or as a branch of "non-accelerator particle physics" , searching for 
the possible imprints of high-energy particle effects in the low-energy realm of gravity. Whether 
current or proposed experiments can actually probe these phenomena meaningfully is an open 
question at the moment, largely because of a dearth of firm theoretical predictions. 

2.3.1 The "fifth" force 

On the phenomenological side, the idea of using EEP tests in this way may have originated in the 
middle 1980s, with the search for a "fifth" force. In 1986, as a result of a detailed reanalysis of 
Eotvos' original data, Fischbach et al. |109j suggested the existence of a fifth force of nature, with 
a strength of about a percent that of gravity, but with a range (as defined by the range A of a 
Yukawa potential, e~ r / A /r) of a few hundred meters. This proposal dovetailed with earlier hints 
of a deviation from the inverse-square law of Newtonian gravitation derived from measurements 
of the gravity profile down deep mines in Australia, and with emerging ideas from particle physics 
suggesting the possible presence of very low-mass particles with gravitational-strength couplings. 
During the next four years numerous experiments looked for evidence of the fifth force by searching 
for composition-dependent differences in acceleration, with variants of the Eotvos experiment or 
with free-fall Galileo- type experiments. Although two early experiments reported positive evidence, 
the others all yielded null results. Over the range between one and 10 4 meters, the null experiments 
produced upper limits on the strength of a postulated fifth force between 10 -3 and 10 -6 of the 
strength of gravity. Interpreted as tests of WEP (corresponding to the limit of infinite-range forces), 
the results of two representative experiments from this period, the free-fall Galileo experiment and 
the early Eot-Wash experiment, are shown in Figure^ At the same time, tests of the inverse-square 
law of gravity were carried out by comparing variations in gravity measurements up tall towers or 
down mines or boreholes with gravity variations predicted using the inverse square law together 
with Earth models and surface gravity data mathematically "continued" up the tower or down the 
hole. Despite early reports of anomalies, independent tower, borehole and seawater measurements 
ultimately showed no evidence of a deviation. Analyses of orbital data from planetary range 
measurements, lunar laser ranging, and laser tracking of the LAGEOS satellite verified the inverse- 
square law to parts in 10 8 over scales of 10 3 to 10 5 km, and to parts in 10 9 over planetary scales 
of several astronomical units |251| . A consensus emerged that there was no credible experimental 
evidence for a fifth force of nature, of a type and range proposed by Fischbach et al. For reviews 
and bibliographies of this episode, see [TH^l HTHl ITTT1 IH1 1275 ] . 

2.3.2 Short-range modifications of Newtonian gravity 

Although the idea of an intermediate-range violation of Newton's gravitational law was dropped, 
new ideas emerged to suggest the possibility that the inverse-square law could be violated at very 
short ranges, below the centimeter range of existing laboratory verifications of the 1/r 2 behavior. 
One set of ideas (151 1131 12221 1221) posited that some of the extra spatial dimensions that come with 
string theory could extend over macroscopic scales, rather than being rolled up at the Planck scale 
of 10" 33 cm, which was then the conventional viewpoint. On laboratory distances large compared 
to the relevant scale of the extra dimension, gravity would fall off as the inverse square, whereas 
on short scales, gravity would fall off as l/i? 2+n , where n is the number of large extra dimensions. 
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Many models favored n — 1 or n = 2. Other possibilities for effective modifications of gravity at 
short range involved the exchange of light scalar particles. 

Following these proposals, many of the high-precision, low-noise methods that were developed 
for tests of WEP were adapted to carry out laboratory tests of the inverse square law of Newto- 
nian gravitation at millimeter scales and below. The challenge of these experiments has been to 
distinguish gravitation-like interactions from electromagnetic and quantum mechanical (Casimir) 
effects. No deviations from the inverse square law have been found to date at distances between 
fO/im and 10 mm [T72*l 11311 [131)1 1541 1171) . For a comprehensive review of both the theory and the 
experiments, see [5]- 

3 Tests of Post-Newtonian Gravity 

3.1 Metric theories of gravity and the strong equivalence principle 

3.1.1 Universal coupling and the metric postulates 

The empirical evidence supporting the Einstein equivalence principle, discussed in the previous 
section, supports the conclusion that the only theories of gravity that have a hope of being viable 
are metric theories, or possibly theories that are metric apart from very weak or short-range 
non- metric couplings (as in string theory). Therefore for the remainder of this review, we shall 
turn our attention exclusively to metric theories of gravity, which assume that (i) there exists a 
symmetric metric, (ii) test bodies follow geodesies of the metric, and (hi) in local Lorentz frames, 
the non-gravitational laws of physics are those of special relativity. 

The property that all non-gravitational fields should couple in the same manner to a single 
gravitational field is sometimes called "universal coupling" . Because of it, one can discuss the 
metric as a property of spacetime itself rather than as a field over spacetime. This is because its 
properties may be measured and studied using a variety of different experimental devices, composed 
of different non-gravitational fields and particles, and, because of universal coupling, the results 
will be independent of the device. Thus, for instance, the proper time between two events is a 
characteristic of spacetime and of the location of the events, not of the clocks used to measure it. 

Consequently, if EEP is valid, the non-gravitational laws of physics may be formulated by taking 
their special relativistic forms in terms of the Minkowski metric r) and simply "going over" to new 
forms in terms of the curved spacetime metric g, using the mathematics of differential geometry. 
The details of this "going over" can be found in standard textbooks ( [1901 1270] . TEGP 3.2. |281| \ 

3.1.2 The strong equivalence principle 

In any metric theory of gravity, matter and non-gravitational fields respond only to the spacetime 
metric g. In principle, however, there could exist other gravitational fields besides the metric, such 
as scalar fields, vector fields, and so on. If, by our strict definition of metric theory, matter does 
not couple to these fields, what can their role in gravitation theory be? Their role must be that 
of mediating the manner in which matter and non-gravitational fields generate gravitational fields 
and produce the metric; once determined, however, the metric alone acts back on the matter in 
the manner prescribed by EEP. 

What distinguishes one metric theory from another, therefore, is the number and kind of 
gravitational fields it contains in addition to the metric, and the equations that determine the 
structure and evolution of these fields. From this viewpoint, one can divide all metric theories of 
gravity into two fundamental classes: "purely dynamical" and "prior-geometric" . 

By "purely dynamical metric theory" we mean any metric theory whose gravitational fields 
have their structure and evolution determined by coupled partial differential field equations. In 
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other words, the behavior of each field is influenced to some extent by a coupling to at least one 
of the other fields in the theory. By "prior geometric" theory we mean any metric theory that 
contains "absolute elements" , fields or equations whose structure and evolution are given a priori, 
and are independent of the structure and evolution of the other fields of the theory. These "absolute 
elements" typically include flat background metrics rj or cosmic time coordinates t. 

General relativity is a purely dynamical theory since it contains only one gravitational field, 
the metric itself, and its structure and evolution are governed by partial differential equations 
(Einstein's equations). Brans-Dicke theory and its generalizations are purely dynamical theories; 
the field equation for the metric involves the scalar field (as well as the matter as source), and 
that for the scalar field involves the metric. Rosen's bimctric theory is a prior-geometric theory: It 
has a non-dynamical, Ricmann-flat background metric r], and the field equations for the physical 
metric g involve r\. 

By discussing metric theories of gravity from this broad point of view, it is possible to draw some 
general conclusions about the nature of gravity in different metric theories, conclusions that are 
reminiscent of the Einstein equivalence principle, but that are subsumed under the name "strong 
equivalence principle" . 

Consider a local, freely falling frame in any metric theory of gravity. Let this frame be small 
enough that inhomogcncities in the external gravitational fields can be neglected throughout its 
volume. On the other hand, let the frame be large enough to encompass a system of gravitating 
matter and its associated gravitational fields. The system could be a star, a black hole, the solar 
system or a Cavendish experiment. Call this frame a "quasi-local Lorentz frame" . To determine 
the behavior of the system we must calculate the metric. The computation proceeds in two stages. 
First we determine the external behavior of the metric and gravitational fields, thereby establishing 
boundary values for the fields generated by the local system, at a boundary of the quasi-local frame 
"far" from the local system. Second, we solve for the fields generated by the local system. But 
because the metric is coupled directly or indirectly to the other fields of the theory, its structure 
and evolution will be influenced by those fields, and in particular by the boundary values taken on 
by those fields far from the local system. This will be true even if we work in a coordinate system in 
which the asymptotic form of in the boundary region between the local system and the external 
world is that of the Minkowski metric. Thus the gravitational environment in which the local 
gravitating system resides can influence the metric generated by the local system via the boundary 
values of the auxiliary fields. Consequently, the results of local gravitational experiments may 
depend on the location and velocity of the frame relative to the external environment. Of course, 
local non-gravitational experiments are unaffected since the gravitational fields they generate are 
assumed to be negligible, and since those experiments couple only to the metric, whose form can 
always be made locally Minkowskian at a given spacetime event. Local gravitational experiments 
might include Cavendish experiments, measurement of the acceleration of massive self-gravitating 
bodies, studies of the structure of stars and planets, or analyses of the periods of "gravitational 
clocks" . We can now make several statements about different kinds of metric theories. 

(i) A theory which contains only the metric g yields local gravitational physics which is inde- 
pendent of the location and velocity of the local system. This follows from the fact that the only 
field coupling the local system to the environment is g, and it is always possible to find a coordinate 
system in which g takes the Minkowski form at the boundary between the local system and the 
external environment (neglecting inhomogeneities in the external gravitational field). Thus the 
asymptotic values of g^ v are constants independent of location, and are asymptotically Lorentz 
invariant, thus independent of velocity. General relativity is an example of such a theory. 

(ii) A theory which contains the metric g and dynamical scalar fields ifA yields local gravi- 
tational physics which may depend on the location of the frame but which is independent of the 
velocity of the frame. This follows from the asymptotic Lorentz invariance of the Minkowski metric 
and of the scalar fields, but now the asymptotic values of the scalar fields may depend on the loca- 
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tion of the frame. An example is Brans-Dickc theory, where the asymptotic scalar field determines 
the effective value of the gravitational constant, which can thus vary as ip varies. On the other 
hand, a form of velocity dependence in local physics can enter indirectly if the asymptotic values 
of the scalar field vary with time cosmologically. Then the rate of variation of the gravitational 
constant could depend on the velocity of the frame. 

(iii) A theory which contains the metric g and additional dynamical vector or tensor fields or 
prior-geometric fields yields local gravitational physics which may have both location and velocity- 
dependent effects. 

These ideas can be summarized in the strong equivalence principle (SEP), which states that: 

1. WEP is valid for self-gravitating bodies as well as for test bodies. 

2. The outcome of any local test experiment is independent of the velocity of the (freely falling) 
apparatus. 

3. The outcome of any local test experiment is independent of where and when in the universe 
it is performed. 

The distinction between SEP and EEP is the inclusion of bodies with self-gravitational interactions 
(planets, stars) and of experiments involving gravitational forces (Cavendish experiments, gravime- 
ter measurements). Note that SEP contains EEP as the special case in which local gravitational 
forces are ignored. 

The above discussion of the coupling of auxiliary fields to local gravitating systems indicates 
that if SEP is strictly valid, there must be one and only one gravitational field in the universe, the 
metric g. These arguments are only suggestive however, and no rigorous proof of this statement is 
available at present. Empirically it has been found that almost every metric theory other than GR 
introduces auxiliary gravitational fields, either dynamical or prior geometric, and thus predicts 
violations of SEP at some level (here we ignore quantum-theory inspired modifications to GR 
involving "i? 2 " terms). The one exception is Nordstrom's 1913 conformally-flat scalar theory |196j . 
which can be written purely in terms of the metric; the theory satisfies SEP, but unfortunately 
violates experiment by predicting no deflection of light. General relativity seems to be the only 
viable metric theory that embodies SEP completely. In Sec. 13.61 we shall discuss experimental 
evidence for the validity of SEP. 

3.2 The parametrized post-Newtonian formalism 

Despite the possible existence of long-range gravitational fields in addition to the metric in var- 
ious metric theories of gravity, the postulates of those theories demand that matter and non- 
gravitational fields be completely oblivious to them. The only gravitational field that enters the 
equations of motion is the metric g. The role of the other fields that a theory may contain can 
only be that of helping to generate the spacetime curvature associated with the metric. Matter 
may create these fields, and they plus the matter may generate the metric, but they cannot act 
back directly on the matter. Matter responds only to the metric. 

Thus the metric and the equations of motion for matter become the primary entities for calcu- 
lating observable effects, and all that distinguishes one metric theory from another is the particular 
way in which matter and possibly other gravitational fields generate the metric. 

The comparison of metric theories of gravity with each other and with experiment becomes 
particularly simple when one takes the slow-motion, weak-field limit. This approximation, known 
as the post- Newtonian limit, is sufficiently accurate to encompass most solar-system tests that 
can be performed in the foreseeable future. It turns out that, in this limit, the spacetime metric 
g predicted by nearly every metric theory of gravity has the same structure. It can be written 
as an expansion about the Minkowski metric (j]^ = diag(— 1, 1, 1, 1)) in terms of dimcnsionless 
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Table 2: The PPN Parameters and their significance (note that has been shown twice to indicate 
that it is a measure of two effects). 



gravitational potentials of varying degrees of smallness. These potentials are constructed from the 
matter variables (Box[5J in imitation of the Newtonian gravitational potential 

U(x,t) = J />(x , ,t)|x-x'|- 1 d?V- (29) 

The "order of smallness" is determined according to the rules U ~ v 2 ~ n ~ p/p ~ e, v % ~ 
\d/dt\/\d/dx\ ~ e 1 / 2 , and so on (we use units in which G = c = 1; see Box 01. 

A consistent post- Newtonian limit requires determination of goo correct through 0(e 2 ), got 
through 0(e 3 / 2 ) and g %j through 0(e) (for details see TEGP 4.1 |2BT] ). The only way that one 
metric theory differs from another is in the numerical values of the coefficients that appear in front 
of the metric potentials. The parametrized post-Newtonian (PPN) formalism inserts parameters 
in place of these coefficients, parameters whose values depend on the theory under study. In the 
current version of the PPN formalism, summarized in Box [21 ten parameters are used, chosen 
in such a manner that they measure or indicate general properties of metric theories of gravity 
(Table |2J) . Under reasonable assumptions about the kinds of potentials that can be present at 
post-Newtonian order (basically only Poisson-like potentials), one finds that ten PPN parameters 
exhaust the possibilities. 

The parameters 7 and (3 are the usual Eddington-Robertson-Schiff parameters used to describe 
the "classical" tests of GR, and are in some sense the most important; they are the only non- 
zero parameters in GR and scalar-tensor gravity. The parameter £ is non-zero in any theory 
of gravity that predicts preferred-location effects such as a galaxy-induced anisotropy in the local 
gravitational constant Gl (also called "Whitehead" effects); a±, ct2, measure whether or not the 
theory predicts post-Newtonian preferred- frame effects; Cl) C2, C3; C4 measure whether or not 
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the theory predicts violations of global conservation laws for total momentum. In TableElwe show 
the values these parameters take (i) in GR, (ii) in any theory of gravity that possesses conservation 
laws for total momentum, called "semi-conservative" (any theory that is based on an invariant 
action principle is semi-conservative), and (iii) in any theory that in addition possesses six global 
conservation laws for angular momentum, called "fully conservative" (such theories automatically 
predict no post-Newtonian preferred- frame effects) . Semi-conservative theories have five free PPN 
parameters (7, f3, £, at, a%) while fully conservative theories have three (7, /3, £). 

The PPN formalism was pioneered by Kenneth Nordtvedt |198) . who studied the post-Newtonian 
metric of a system of gravitating point masses, extending earlier work by Eddington, Robertson 
and Schiff (TEGP 4.2 |2HU)- Will [23] generalized the framework to perfect fluids. A general 
and unified version of the PPN formalism was developed by Will and Nordtvedt. The canonical 
version, with conventions altered to be more in accord with standard textbooks such as |19()j . is 
discussed in detail in TEGP 4 |281) . Other versions of the PPN formalism have been developed 
to deal with point masses with charge, fluid with anisotropic stresses, bodies with strong internal 
gravity, and post-post-Newtonian effects (TEGP 4.2, 14.2 |2"5T]). 
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Box 2. The Parametrized Post-Newtonian Formalism 



1. Coordinate system: The framework uses a nearly globally Lorentz coordinate system 
in which the coordinates are (t . X « X f X ). Three-dimensional, Euclidean vector notation 
is used throughout. All coordinate arbitrariness ("gauge freedom") has been removed by 
specialization of the coordinates to the standard PPN gauge (TEGP 4.2 |281j ). Units are 
chosen so that G = c = 1, where G is the physically measured Newtonian constant far from 
the solar system. 

2. Matter variables: 

• p: density of rest mass as measured in a local freely falling frame momentarily comoving 
with the gravitating matter; 

• v l = (dx l /dt): coordinate velocity of the matter; 

• w l : coordinate velocity of the PPN coordinate system relative to the mean rest-frame 
of the universe; 

• p: pressure as measured in a local freely falling frame momentarily comoving with the 
matter; 

• II: internal energy per unit rest mass. It includes all forms of non-rest-mass, non- 
gravitational energy, e.g. energy of compression and thermal energy. 

3. PPN parameters: 7, @, f, a x , a 2 , ar 3 , C2, C3, Ci- 

4. Metric: 

3oo = -l + 2U- 2/3C/ 2 - 2& w + (2 7 + 2 + a 3 + Ci ~ 2£)*i 

+2(3 7 - 2[3 + 1 + C2 + 0*2 + 2(1 + C 3 )$3 + 2(3 7 + 3C 4 - 2£)$ 4 
-(Ci - 2£)A - (ai - a 2 - a 3 )w 2 U - a 2 w l w 3 U lJ + (2a 3 - ai)w l Vi 
+0(e 3 ), 

goi = ~{4rf + 3 + ai - a 2 + Ci - 2£)V^ - |(l + ar 2 -Ci + 2£)TT, 

~-(a x - 2a 2 )w l U - a 2 wW l0 + 0(e 5/2 ), 
9ij = (1 + 2jU)Sij + 0(e 2 ). 

5. Metric potentials: 

U 



<f> w 
A 
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p'(x-x')i(x-x') 3 ; 3 , 

|x-x'| 3 

p/(x-x') /x'-x" x-x" 



p'[v' ■ (x - x')] 1 



'„/2 



p V 



,/|3 



d 3 x'd 3 x" 



J |x-x'| 
7 |x-x'| 

7 |x-x'| 3 

6. Stress— energy tensor (perfect fluid): 

T 00 = p(l + Il + v 2 + 2U), 

T 0i = pv i (l + U + v 2 + 2U + p/p), 

T i] = pt;V(l + n + v 2 + 2U + p/p) +p5 ij (l - 2-yU). 

7. Equations of motion: 

• Stressed matter: T» v . v = 0, 

• Test bodies: d 2 x^/d\ 2 + T^ uX (dx v /d\)(dx x /d\) = 0, 

• Maxwell's equations: F^" . v = inJ^, F^ = A v .^ - A^ v . 
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3.3 Competing theories of gravity 



One of the important applications of the PPN formalism is the comparison and classification of 
alternative metric theories of gravity. The population of viable theories has fluctuated over the 
years as new effects and tests have been discovered, largely through the use of the PPN framework, 
which eliminated many theories thought previously to be viable. The theory population has also 
fluctuated as new, potentially viable theories have been invented. 

In this review, we shall focus on GR, the general class of scalar-tensor modifications of it, of 
which the Jordan-Fierz-Brans-Dicke theory (Brans-Dicke, for short) is the classic example, and 
vector-tensor theories. The reasons are several-fold: 

• A full compendium of alternative theories circa 1981 is given in TEGP 5 |281| . 

• Many alternative metric theories developed during the 1970s and 1980s could be viewed as 
"straw-man" theories, invented to prove that such theories exist or to illustrate particular 
properties. Few of these could be regarded as well-motivated theories from the point of view, 
say, of field theory or particle physics. 

• A number of theories fall into the class of "prior-geometric" theories, with absolute elements 
such as a flat background metric in addition to the physical metric. Most of these theories 
predict "preferred-frame" effects, that have been tightly constrained by observations (see 
Sec. I3.6.2")l . An example is Rosen's bimetric theory. 

• A large number of alternative theories of gravity predict gravitational wave emission sub- 
stantially different from that of general relativity, in strong disagreement with observations 
of the binary pulsar (see Sec. EJ. 

• Scalar-tensor modifications of GR have become very popular in unification schemes such as 
string theory, and in cosmological model building. Because the scalar fields could be massive, 
the potentials in the post-Newtonian limit could be modified by Yukawa-like terms. 

• Vector-tensor theories have attracted recent attention, in the spirit of the SME fSec. I2.2!lj) . 
as models for violations of Lorentz invariance in the gravitational sector. 

3.3.1 General relativity 

The metric g is the sole dynamical field, and the theory contains no arbitrary functions or parame- 
ters, apart from the value of the Newtonian coupling constant G, which is measurable in laboratory 
experiments. Throughout this article, we ignore the cosmological constant Ac- We do this despite 
recent evidence, from supernova data, of an accelerating universe, which would indicate either a 
non-zero cosmological constant or a dynamical "dark energy" contributing about 70 percent of 
the critical density. Although Ac has significance for quantum field theory, quantum gravity, and 
cosmology, on the scale of the solar-system or of stellar systems its effects are negligible, for the 
values of Ac inferred from supernova observations. 

The field equations of GR are derivable from an invariant action principle 51 = 0, where 



where R is the Ricci scalar, and I m is the matter action, which depends on matter fields ip m 
universally coupled to the metric g. By varying the action with respect to t^,„, we obtain the field 
equations 




(30) 



G>„ = R^v — -^g^vR — 8ttGT^ 



(31) 
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where T^ v is the matter energy-momentum tensor. General covariance of the matter action implies 
the equations of motion - v = 0; varying 7 m with respect to V'm yields the matter field equations 
of the Standard Model. By virtue of the absence of prior-geometric elements, the equations of 
motion are also a consequence of the field equations via the Bianchi identities ' = 0. 

The general procedure for deriving the post-Newtonian limit of metric theories is spelled out in 
TEGP 5.1 [2HD, and is described in detail for GR in TEGP 5.2 [2HD- The PPN parameter values 
are listed in Table 

3.3.2 Scalar-tensor theories 

These theories contain the metric g, a scalar held tp, a potential function V(ip), and a coupling 
function A(<p) (generalizations to more than one scalar field have also been carried out |76j). 
For some purposes, the action is conveniently written in a non-metric representation, sometimes 
denoted the "Einstein frame" , in which the gravitational action looks exactly like that of GR: 

I = (mirG)- 1 f[R- Tg^d^dvV - V(<p)](-g) 1/2 d 4 x + I m (^ m ,A 2 (p)g^), (32) 

where R = g^ v i? M „ is the Ricci scalar of the "Einstein" metric g^. (Apart from the scalar potential 
term V(p), this corresponds to Eq. with G(cp) = (AirGy 1 , U(ip) = 1, and M(tp) oc A(cp).) 
This representation is a "non-metric" one because the matter fields ip m couple to a combination 
of ip and g^v Despite appearances, however, it is a metric theory, because it can be put into a 
metric representation by identifying the "physical metric" 

g^ ee A 2 (if)g^. (33) 

The action can then be rewritten in the metric form 

/ = (IGttG)- 1 f [<t>R - (TM^y^cU - tfVK-gf/Wx + I m Wm,9» v ), (34) 
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where 



ct> = A^)-\ 

3 + 2co(0) = a(^)- 2 , (35) 
a(<p) = d(kiA((p))/dip. 

The Einstein frame is useful for discussing general characteristics of such theories, and for some 
cosmological applications, while the metric representation is most useful for calculating observable 
effects. The field equations, post-Newtonian limit and PPN parameters are discussed in TEGP 
5.3 |281j . and the values of the PPN parameters are listed in Table 
The parameters that enter the post-Newtonian limit are 

u = u(<M, A=[(dcj/d<P)(3 + 2Lu)- 2 (4 + 2uj)-% , (36) 

where 4>o is the value of <f> today far from the system being studied, as determined by appropriate 
cosmological boundary conditions. In Brans-Dicke theory (u{<p) = o^bd = constant), the larger 
the value of wbd, the smaller the effects of the scalar field, and in the limit cjbd — ► °o (ckq — * 0), 
the theory becomes indistinguishable from GR in all its predictions. In more general theories, the 
function u>{4>) could have the property that, at the present epoch, and in weak-field situations, the 
value of the scalar field <po is such that u> is very large and A is very small (theory almost identical 
to GR today) , but that for past or future values of <fi, or in strong-field regions such as the interiors 
of neutron stars, u> and A could take on values that would lead to significant differences from GR. 
It is useful to point out that all versions of scalar-tensor gravity predict that 7 < 1 (Table 0. 

Damour and Esposito-Farese |76| have adopted an alternative parametrization of scalar-tensor 
theories, in which one expands hiA(ip) about a cosmological background field value ipg: 

In A(<p) = a ((p - <Po) + xPo(<P ~ fo f + ■ • ■ ■ (37) 

A precisely linear coupling function produces Brans-Dicke theory, with — 1/(2cjbd + 3), or 
1/(2 + cjbd) = 2o!q/ (1 + q.q). The function In A(ip) acts as a potential for the scalar field ip within 
matter, and, if /3q > 0, then during cosmological evolution, the scalar field naturally evolves toward 
the minimum of the potential, i.e. toward ao ~ 0, lu — > 00, or toward a theory close to, though not 
precisely GR [831 184) . Estimates of the expected relic deviations from GR today in such theories 
depend on the cosmological model, but range from 1CP 5 to a few times 10 -7 for I7 — 1|. 

Negative values of (3q correspond to a "locally unstable" scalar potential (the overall theory is 
still stable in the sense of having no tachyons or ghosts). In this case, objects such as neutron stars 
can experience a "spontaneous scalarization" , whereby the interior values of tp can take on values 
very different from the exterior values, through non-linear interactions between strong gravity and 
the scalar field, dramatically affecting the stars' internal structure and leading to strong violations 
of SEP. On the other hand, in the case (3$ < 0, one must confront that fact that, with an unstable 
ip potential, cosmological evolution would presumably drive the system away from the peak where 
ao w 0, toward parameter values that could be excluded by solar system experiments. 

Scalar fields coupled to gravity or matter are also ubiquitous in particle-physics-inspired models 
of unification, such as string theory [2551 11771 1881 1851 I86| . In some models, the coupling to matter 
may lead to violations of EEP, which could be tested or bounded by the experiments described 
in Sec. 12.11 In many models the scalar field could be massive; if the Compton wavelength is of 
macroscopic scale, its effects are those of a "fifth force" . Only if the theory can be cast as a metric 
theory with a scalar field of infinite range or of range long compared to the scale of the system in 
question (solar system) can the PPN framework be strictly applied. If the mass of the scalar field 
is sufficiently large that its range is microscopic, then, on solar-system scales, the scalar field is 
suppressed, and the theory is essentially equivalent to general relativity. 
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3.3.3 Vector- Tensor theories 



These theories contain the metric g and a dynamical, typically timelike, four- vector field u^ L . In 
some models, the four-vector is unconstrained, while in others, called Einstein- iEther theories it 
is constrained to be timelike with unit norm. The most general action for such theories that is 
quadratic in derivatives of the vector is given by 

(-g) 1/2 d 4 2; + / m O m ,5 Mt ,), , 

(38) 
where 

Kp = cig^g a p + c 2 S^Sf 3 + c 3 6%6» - c^u"g a p . (39) 

The coefficients Cj are arbitrary. In the unconstrained theories, A = and lj is arbitrary. In the 
constrained theories, A is a Lagrange multiplier, and by virtue of the constraint u^u 11 — — 1, the 
factor uju^u^ in front of the Ricci scalar can be absorbed into a rescaling of G; equivalently, in the 
constrained theories, we can set w = 0. Note that the possible term u^u^R^ can be shown under 
integration by parts to be equivalent to a linear combination of the terms involving c 2 and C3 . 

Unconstrained theories were studied during the 1970s as "straw-man" alternatives to GR. In 
addition to having up to four arbitrary parameters, they also left the magnitude of the vector field 
arbitrary, since it satisfies a linear homogenous vacuum field equation of the form C,u^ = (04 = 
in all such cases studied). Indeed, this latter fact was one of most serious defects of these theories. 
Each theory studied corresponds to a special case of the action l)38[l . all with A = 0: 

• General vector-tensor theory (u, t, e, rj), TEGP 5.4 |281| : The gravitational Lagrangian for 
this class of theories had the form R + uiu^u^R + rju^u^R^ — eF^F^ +tV m Ui / V a 'u 1/ , where 
Ffj, u = Vfj,u u — V„w M , corresponding to the values c\ = 2e — t, c 2 = — 77, C\ + c 2 + c 3 = — r, 
C4 = 0. In these theories 7, (3, a,\ and a 2 are complicated functions of the parameters and of 
u 2 = —u^u^, while the rest vanish. 

• Will-Nordtvedt theory |290| : This is the special case ci = — 1, C2 = C3 = C4 = 0. In this 
theory, the PPN parameters are given by 7 = f3 = 1, «2 = u 2 /(l + u 2 /2), and zero for the 
rest. 

• Hellings-Nordtvedt theory (oj) |129) : This is the special case C\ — 2, c 2 = 2oj, ci + c 2 + c 3 = 
= c 4 . Here 7, (3, ct\ and a 2 are complicated functions of the parameters and of it 2 , while 
the rest vanish. 

The Einstein-iEther theories were motivated in part by a desire to explore possibilities for 
violations of Lorentz invariance in gravity, in parallel with similar studies in matter interactions, 
such as the SME. The general class of theories was analysed by Jacobson and collaborators [1381 
HHH EM IMl EH , motivated in part by \T57\ . 

Analysing the post-Newtonian limit, they were able to infer values of the PPN parameters 7 
and P as follows |114| : 



/ = (16ttG)- 



(1+luu^)R- 



1) 



7=1, 0=1, 
£ = a> 3 = Ci = C2 = C3 = Ci , 
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ai = - 
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(2ci 3 - cu) 2 
ci 23 (2 - cu) 
12c 3 ci 3 + 2c 2 ci 4 (l 



2cu) + (c? - c§)(4 - 6c 13 + 7c M ) 



(2-c u )(2ci -c\ + cl) 



(43) 



where C123 = c\ + c 2 + C3, C13 = c\ — C3, C14 = Ci — C4, subject to the constraints ci2 3 7^ 0, C14 7^ 2, 



2ci 



7^0. 



By requiring that gravitational- wave modes have real (as opposed to imaginary) frequencies, 
one can impose the bounds c\j(c\ + C4) > and (ci + c 2 + C3)/(ci + C4) > 0. Considerations of 
positivity of energy impose the constraints c\ > 0, c\ + C4 > and c\ + c 2 + c 3 > 0. 

3.4 Tests of the parameter 7 

With the PPN formalism in hand, we are now ready to confront gravitation theories with the 
results of solar-system experiments. In this section we focus on tests of the parameter 7, consisting 
of the deflection of light and the time delay of light. 

3.4.1 The deflection of light 

A light ray (or photon) which passes the Sun at a distance d is deflected by an angle 

1 



= -(l+7)(4m /d)[(l + coB*)/2] 



(44) 



(TEGP 7.1 |281] V where is the mass of the Sun and $ is the angle between the Earth-Sun line 
and the incoming direction of the photon (Figure For a grazing ray, d ~ <i©, $ « 0, and 



-(1+ 7 )1."7505, 



(45) 



independent of the frequency of light. Another, more useful expression gives the change in the 
relative angular separation between an observed source of light and a nearby reference source as 
both rays pass near the Sun: 



;(l + 7) 



4m r . 



• cos x 



477J G 



1 + cos $ r 



(46) 



where d and d r are the distances of closest approach of the source and reference rays respectively, 
$ r is the angular separation between the Sun and the reference source, and x is the angle between 
the Sun-source and the Sun- reference directions, projected on the plane of the sky (Figure 0J. 
Thus, for example, the relative angular separation between the two sources may vary if the line of 
sight of one of them passes near the Sun (d ~ Rq, d r ^> d, x varying with time). 

It is interesting to note that the classic derivations of the deflection of light that use only the 
corpuscular theory of light (Cavendish 1784, von Soldner 1803 |277| 1. or the principle of equivalence 
(Einstein 1911), yield only the "1/2" part of the coefficient in front of the expression in Eq. I|44|l . 
But the result of these calculations is the deflection of light relative to local straight lines, as estab- 
lished for example by rigid rods; however, because of space curvature around the Sun, determined 
by the PPN parameter 7, local straight lines are bent relative to asymptotic straight lines far from 
the Sun by just enough to yield the remaining factor "7/2". The first factor "1/2" holds in any 
metric theory, the second "7/2" varies from theory to theory. Thus, calculations that purport to 
derive the full deflection using the equivalence principle alone are incorrect. 

The prediction of the full bending of light by the Sun was one of the great successes of Einstein's 
GR. Eddington's confirmation of the bending of optical starlight observed during a solar eclipse in 



33 



Figure 4: Geometry of light deflection measurements. 



the first days following World War 1 helped make Einstein famous. However, the experiments of 
Eddington and his co-workers had only 30 percent accuracy, and succeeding experiments were not 
much better: The results were scattered between one half and twice the Einstein value (Figure^, 
and the accuracies were low. 

However, the development of radio-interferometery, and later of VLBI, very-long-baseline ra- 
dio interferometry, produced greatly improved determinations of the deflection of light. These 
techniques now have the capability of measuring angular separations and changes in angles to 
accuracies better than 100 microarcseconds. Early measurements took advantage of a series of 
heavenly coincidences: Each year, groups of strong quasistellar radio sources pass very close to 
the Sun (as seen from the Earth), including the group 3C273, 3C279, and 3C48, and the group 
0111+02, 0119+11 and 0116+08. As the Earth moves in its orbit, changing the lines of sight of the 
quasars relative to the Sun, the angular separation 69 between pairs of quasars varies (Eq. I|46[) ). 
The time variation in the quantities d, d r , \ and <I> r in Eq. I|46[) is determined using an accurate 
ephemcris for the Earth and initial directions for the quasars, and the resulting prediction for 86 
as a function of time is used as a basis for a least-squares fit of the measured 59, with one of 
the fitted parameters being the coefficient i(l + 7). A number of measurements of this kind over 
the period 1969-1975 yielded an accurate determination of the coefficient 4(1+7). A 1995 VLBI 
measurement using 3C273 and 3C279 yielded (l + 7)/2 = 0.9996+0.0017 |TS5] . 

In recent years, transcontinental and intercontinental VLBI observations of quasars and radio 
galaxies have been made primarily to monitor the Earth's rotation ("VLBI" in Figure [SJ. These 
measurements are sensitive to the deflection of light over almost the entire celestial sphere (at 90° 
from the Sun, the deflection is still 4 milliarcseconds) . A 2004 analysis of almost 2 million VLBI 
observations of 541 radio sources, made by 87 VLBI sites yielded (1 + 7 )/2 = 0.99992 ± 0.00023, 
or equivalently, 7 - 1 = (-1.7 ± 4.5) x 10" 4 

Analysis of observations made by the Hipparcos optical astrometry satellite yielded a test at 
the level of 0.3 percent |116) . A VLBI measurement of the deflection of light by Jupiter was 
reported; the predicted deflection of about 300 microarcseconds was seen with about 50 percent 
accuracy 258 . The results of light-deflection measurements are summarized in Figure [S] 
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Figure 5: Measurements of the coefficient (1 + j)/2 from light deflection and time delay measure- 
ments. Its GR value is unity. The arrows at the top denote anomalously large values from early 
eclipse expeditions. The Shapiro time-delay measurements using the Cassini spacecraft yielded an 
agreement with GR to 10 -3 percent, and VLBI light deflection measurements have reached 0.02 
percent. Hipparcos denotes the optical astrometry satellite, which reached 0.1 percent. 
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3.4.2 The time delay of light 

A radar signal sent across the solar system past the Sun to a planet or satellite and returned to 
the Earth suffers an additional non-Newtonian delay in its round-trip travel time, given by (see 
Figure El 

St = 2(1 + 7 )m ln[(r e + x ffi • n)(r e - x c • n)/d 2 ], (47) 

where x e (x®) are the vectors, and r c (r ffi ) are the distances from the Sun to the source (Earth), 
respectively (TEGP 7.2 [281] ). For a ray which passes close to the Sun, 

6t& -(1 + 7 )[240- 201n(d 2 /r)] (is, (48) 

where d is the distance of closest approach of the ray in solar radii, and r is the distance of the 
planet or satellite from the Sun, in astronomical units. 

In the two decades following Irwin Shapiro's 1964 discovery of this effect as a theoretical con- 
sequence of GR, several high-precision measurements were made using radar ranging to targets 
passing through superior conjunction. Since one does not have access to a "Newtonian" signal 
against which to compare the round-trip travel time of the observed signal, it is necessary to do a 
differential measurement of the variations in round-trip travel times as the target passes through 
superior conjunction, and to look for the logarithmic behavior of Eq. I|48[l. In order to do this 
accurately however, one must take into account the variations in round-trip travel time due to the 
orbital motion of the target relative to the Earth. This is done by using radar-ranging (and possibly 
other) data on the target taken when it is far from superior conjunction (i.e. when the time-delay 
term is negligible) to determine an accurate ephemeris for the target, using the ephemeris to pre- 
dict the PPN coordinate trajectory x e (t) near superior conjunction, then combining that trajectory 
with the trajectory of the Earth x s (t) to determine the Newtonian round-trip time and the loga- 
rithmic term in Eq. (|48[) . The resulting predicted round-trip travel times in terms of the unknown 
coefficient |(1 +7) are then fit to the measured travel times using the method of least-squares, 
and an estimate obtained for |(1 + 7). 

The targets employed included planets, such as Mercury or Venus, used as passive reflectors 
of the radar signals ("passive radar"), and artificial satellites, such as Mariners 6 and 7, Voy- 
ager 2, the Viking Mars landers and orbiters, and the Cassini spacecraft to Saturn, used as active 
retransmitters of the radar signals ("active radar"). 

The results for the coefficient ^(1 + 7) of all radar time-delay measurements performed to date 
(including a measurement of the one-way time delay of signals from the millisecond pulsar PSR 
1937+21) are shown in Figure El (see TEGP 7.2 [2HU for discussion and references). The 1976 
Viking experiment resulted in a 0.1 percent measurement 223 . 

A significant improvement was reported in 2003 from Doppler tracking of the Cassini spacecraft 
while it was on its way to Saturn with a result 7 — 1 = (2.1 ± 2.3) x 10~ 5 . This was made 
possible by the ability to do Doppler measurements using both X-band (7175 MHz) and Ka-band 
(34316 MHz) radar, thereby significantly reducing the dispersive effects of the solar corona. In 
addition, the 2002 superior conjunction of Cassini was particularly favorable: with the spacecraft 
at 8.43 astronomical units from the Sun, the distance of closest approach of the radar signals to 
the Sun was only 1.6 Rq. 

From the results of the Cassini experiment, we can conclude that the coefficient i(l + 7) must 
be within at most 0.0012 percent of unity. Scalar-tensor theories must have u> > 40000 to be 
compatible with this constraint. 

3.4.3 Shapiro time delay and the speed of gravity 

In 2001, Kopeikin [152] suggested that a measurement of the time delay of light from a quasar 
as the light passed by the planet Jupiter could be used to measure the speed of the gravitational 
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Table 4: Current limits on the PPN parameters. Here ?7n is a combination of other parameters 
given by m = A3 - 7 - 3 - 10£/3 - ol\ + 2a 2 /3 - 2&/3 - C 2 /3. 



interaction. He argued that, since Jupiter is moving relative to the solar system, and since gravity 
propagates with a finite speed, the gravitational field experienced by the light ray should be affected 
by gravity's speed, since the field experienced at one time depends on the location of the source a 
short time earlier, depending on how fast gravity propagates. According to his calculations, there 
should be a post ^-Newtonian correction to the normal Shapiro time-delay formula (|47|) which 
depends on the velocity of Jupiter and on the velocity of gravity. On September 8, 2002, Jupiter 
passed almost in front of a quasar, and Kopeikin and Fomalont made precise measurements of the 
Shapiro delay with picosecond timing accuracy, and claimed to have measured the correction term 
to about 20 percent [H51 IH8l IT531 H53| . 

However, several authors pointed out that this 1.5PN effect does not depend on the speed 
of propagation of gravity, but rather only depends on the speed of light |16l 12881 12331 1511 1234) . 
Intuitively, if one is working to only first order in v /c, then all that counts is the uniform motion of 
the planet, Jupiter (its acceleration about the Sun contributes a higher-order, unmeasurably small 
effect). But if that is the case, then the principle of relativity says that one can view things from 
the rest frame of Jupiter. In this frame, Jupiter's gravitational field is static, and the speed of 
propagation of gravity is irrelevant. A detailed post-Newtonian calculation of the effect was done 
using a variant of the PPN framework, in a class of theories in which the speed of gravity could be 
different from that of light 288 , and found explicitly that, at first order in v/c, the effect depends 
on the speed of light, not the speed of gravity, in line with intuition. Effects dependent upon the 
speed of gravity show up only at higher order in v/c. Kopeikin gave a number of arguments in 
opposition to this interpretation [1541 1151)1 11491 H5T) . On the other hand, the v/c correction term 
does show a dependence on the PPN parameter a l5 which could be non-zero in theories of gravity 
with a differing speed c g of gravity (see Eq. (7) of [288 ) . But existing tight bounds on ct\ from 
other experiments (Table already far exceed the capability of the Jupiter VLBI experiment. 

3.5 The perihelion shift of Mercury 

The explanation of the anomalous perihelion shift of Mercury's orbit was another of the triumphs 
of GR. This had been an unsolved problem in celestial mechanics for over half a century, since 
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the announcement by Le Verrier in 1859 that, after the perturbing effects of the planets on Mer- 
cury's orbit had been accounted for, and after the effect of the precession of the equinoxes on the 
astronomical coordinate system had been subtracted, there remained in the data an unexplained 
advance in the perihelion of Mercury. The modern value for this discrepancy is 43 arcseconds 
per century. A number of ad hoc proposals were made in an attempt to account for this excess, 
including, among others, the existence of a new planet Vulcan near the Sun, a ring of planetoids, 
a solar quadrupole moment and a deviation from the inverse-square law of gravitation, but none 
was successful. General relativity accounted for the anomalous shift in a natural way without 
disturbing the agreement with other planetary observations. 

The predicted advance per orbit Ad), including both relativistic PPN contributions and the 
Newtonian contribution resulting from a possible solar quadrupole moment, is given by 



where m = mi + m 2 and fj, = mim^/m are the total mass and reduced mass of the two-body 
system respectively; p = a(l — e 2 ) is the semi-latus rectum of the orbit, with the semi-major axis a 
and the eccentricity e; R is the mean radius of the oblate body; and J 2 is a dimensionless measure 
of its quadrupole moment, given by J2 — (C ~ A)/miR 2 , where C and A are the moments of 
inertia about the body's rotation and equatorial axes, respectively (for details of the derivation see 
TEGP 7.3 281 ]). We have ignored preferred-frame and galaxy-induced contributions to Ad); these 
are discussed in TEGP 8.3 |2BT] . 

The first term in Eq. H49[l is the classical relativistic perihelion shift, which depends upon the 
PPN parameters 7 and j3. The second term depends upon the ratio of the masses of the two 
bodies; it is zero in any fully conservative theory of gravity {ct\ = a 2 = 0:3 = C2 = 0); it is also 
negligible for Mercury, since p/m w rriMcrc/m© ~ 2 x 1CP 7 . We shall drop this term henceforth. 

The third term depends upon the solar quadrupole moment J 2 . F° r a Sun that rotates uni- 
formly with its observed surface angular velocity, so that the quadrupole moment is produced by 
centrifugal flattening, one may estimate J 2 to be ~ 1 x 10~ 7 . This actually agrees reasonably well 
with values inferred from rotating solar models that are in accord with observations of the nor- 
mal modes of solar oscillations (helioseismology) ; the latest inversions of helioseismology data give 
J 2 = (2.2 ± 0.1) x 1CT 7 [2IElEniEIIllIEH]- Substituting standard orbital elements and physical 
constants for Mercury and the Sun we obtain the rate of perihelion shift d), in seconds of arc per 



Now, the measured perihelion shift of Mercury is known accurately: after the perturbing effects 
of the other planets have been accounted for, the excess shift is known to about 0.1 percent from 
radar observations of Mercury between 1966 and 1990 239 . Analysis of data taken since 1990 
could improve the accuracy. The solar oblateness effect is smaller than the observational error, so 
we obtain the PPN bound \2*y - j3 - 1| < 3 x 10" 3 . 

3.6 Tests of the strong equivalence principle 

The next class of solar-system experiments that test relativistic gravitational effects may be called 
tests of the strong equivalence principle (SEP). In Sec. 13.1.21 we pointed out that many metric 
theories of gravity (perhaps all except GR) can be expected to violate one or more aspects of 
SEP. Among the testable violations of SEP are a violation of the weak equivalence principle for 
gravitating bodies that leads to perturbations in the Earth-Moon orbit; preferred-location and 



Ad) = 




century, 




(50) 
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preferred-frame effects in the locally measured gravitational constant that could produce observ- 
able geophysical effects; and possible variations in the gravitational constant over cosmological 
timescales. 



3.6.1 The Nordtvedt effect and the lunar Eotvos experiment 

In a pioneering calculation using his early form of the PPN formalism, Nordtvedt |197| showed that 
many metric theories of gravity predict that massive bodies violate the weak equivalence principle 
- that is, fall with different accelerations depending on their gravitational self-energy. Dicke 229 
argued that such an effect would occur in theories with a spatially varying gravitational constant, 
such as scalar-tensor gravity. For a spherically symmetric body, the acceleration from rest in an 
external gravitational potential U has the form 

a = (m p /m)V£7, 

TOp/m = 1 - r^(E s /m), (51) 

10 2 2 1 

??n = 4/? - 7 - 3 - y£ - ai + -a 2 - -Ci - gC2, 

where E g is the negative of the gravitational self-energy of the body (E g > 0). This violation of the 
massive-body equivalence principle is known as the "Nordtvedt effect" . The effect is absent in GR 
(r/N = 0) but present in scalar-tensor theory (77N = 1/(2 +w) + 4A) . The existence of the Nordtvedt 
effect does not violate the results of laboratory Eotvos experiments, since for laboratory-sized 
objects, E g /m < 10^ 27 , far below the sensitivity of current or future experiments. However, for 
astronomical bodies, E g /m may be significant (3.6 x 10~ 6 for the Sun, 10~ 8 for Jupiter, 4.6 x 1CP 10 
for the Earth, 0.2 x 10~ 10 for the Moon). If the Nordtvedt effect is present (??n 7^ 0) then the Earth 
should fall toward the Sun with a slightly different acceleration than the Moon. This perturbation 
in the Earth-Moon orbit leads to a polarization of the orbit that is directed toward the Sun as 
it moves around the Earth-Moon system, as seen from Earth. This polarization represents a 
perturbation in the Earth-Moon distance of the form 

Sr = 13.1 ?7n cos(tt>o — u s )t [m], (52) 

where and w s are the angular frequencies of the orbits of the Moon and Sun around the Earth (see 
TEGP 8.1 |281) for detailed derivations and references; for improved calculations of the numerical 
coefficient, see [203 EU)- 

Since August 1969, when the first successful acquisition was made of a laser signal reflected 
from the Apollo 11 retroreflector on the Moon, the lunar laser-ranging experiment (LLR) has 
made regular measurements of the round-trip travel times of laser pulses between a network of 
observatories and the lunar retroreflectors, with accuracies that are at the 50 ps (1 cm) level, and 
that may soon approach 5 ps (1 mm). These measurements are fit using the method of least-squares 
to a theoretical model for the lunar motion that takes into account perturbations due to the Sun 
and the other planets, tidal interactions, and post-Newtonian gravitational effects. The predicted 
round-trip travel times between retroreflector and telescope also take into account the librations of 
the Moon, the orientation of the Earth, the location of the observatories, and atmospheric effects 
on the signal propagation. The "Nordtvedt" parameter 77N along with several other important 
parameters of the model are then estimated in the least-squares method. 

Numerous ongoing analyses of the data find no evidence, within experimental uncertainty, for 
the Nordtvedt effect [295 296 (for earlier results see |98M294l|l93| ). These results represent a limit 
on a possible violation of WEP for massive bodies of about 1.4 parts in 10 13 (compare Figure 

However, at this level of precision, one cannot regard the results of lunar laser ranging as a 
"clean" test of SEP until one eliminates the possibility of a compensating violation of WEP for 
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the two bodies, because the chemical compositions of the Earth and Moon differ in the relative 
fractions of iron and silicates. To this end, the Eot-Wash group carried out an improved test of 
WEP for laboratory bodies whose chemical compositions mimic that of the Earth and Moon. The 
resulting bound of 1.4 parts in 10 13 |21l 14") from composition effects reduces the ambiguity in the 
lunar laser ranging bound, and establishes the firm SEP test at the level of about 2 parts in 10 13 . 
These results can be summarized by the Nordtvedt parameter bound \rj^\ = (4.4 ± 4.5) x 10~ 4 . 

In the future, the Apache Point Observatory for Lunar Laser-ranging Operation (APOLLO) 
project, a joint effort by researchers from the Universities of Washington, Seattle, and California, 
San Diego, plans to use enhanced laser and telescope technology, together with a good, high- 
altitude site in New Mexico, to improve the lunar laser-ranging bound by as much as an order of 
mag nitude 

In GR, the Nordtvedt effect vanishes; at the level of several centimeters and below, a number 
of non-null general relativistic effects should be present |202) . 

Tests of the Nordtvedt effect for neutron stars have also been carried out using a class of 
systems known as wide-orbit binary millisecond pulsars (WBMSP), which are pulsar- white-dwarf 
binary systems with small orbital eccentricities. In the gravitational field of the galaxy, a non- 
zero Nordtvedt effect can induce an apparent anomalous eccentricity pointed toward the galactic 
center [85], which can be bounded using statistical methods, given enough WBMSPs (see |245) 
for a review and references). Using data from 21 WBMSPs, including recently discovered highly 
circular systems, Stairs et al. |244| obtained the bound A < 5.6 x 10~ 3 , where A = rj-^(E g /M)-^s- 
Because (E g /M)ns ~ 0.1 for typical neutron stars, this bound does not compete with the bound 
on ?7n from Lunar laser ranging; on the other hand, it does test SEP in the strong-field regime 
because of the presence of the neutron stars. 

3.6.2 Preferred-frame and preferred-location effects 

Some theories of gravity violate SEP by predicting that the outcomes of local gravitational experi- 
ments may depend on the velocity of the laboratory relative to the mean rest frame of the universe 
(preferred-frame effects) or on the location of the laboratory relative to a nearby gravitating body 
(preferred-location effects). In the post-Newtonian limit, preferred- frame effects are governed by 
the values of the PPN parameters a\, a>2, and a^, and some preferred- location effects are governed 
by £ (see Table EJ. 

The most important such effects are variations and anisotropies in the locally-measured value 
of the gravitational constant, which lead to anomalous Earth tides and variations in the Earth's 
rotation rate; anomalous contributions to the orbital dynamics of planets and the Moon; self- 
accelerations of pulsars, and anomalous torques on the Sun that would cause its spin axis to be 
randomly oriented relative to the ecliptic (see TEGP 8.2, 8.3, 9.3 and 14.3 (c) |281| 1. An bound 
on (X3 of 4 x 10~ 20 from the period derivatives of 21 millisecond pulsars was reported in |23|544 ; 
improved bounds on a% were achieved using lunar laser ranging data |192| . and using observations 
of the circular binary orbit of the pulsar J2317+1439 [22] ■ Negative searches for these effects have 
produced strong constraints on the PPN parameters (Table HJ. 

3.6.3 Constancy of the Newtonian gravitational constant 

Most theories of gravity that violate SEP predict that the locally measured Newtonian gravitational 
constant may vary with time as the universe evolves. For the scalar-tensor theories listed in TableOl 
the predictions for G/G can be written in terms of time derivatives of the asymptotic scalar field. 
Where G does change with cosmic evolution, its rate of variation should be of the order of the 
expansion rate of the universe, i.e. G/G ~ Ho, where Hq is the Hubble expansion parameter and 
is given by Ho — 100 h km s _1 Mpc -1 = 1.02 x 10~ 10 h yr _1 , where current observations of the 
expansion of the universe give h w 0.73 ± 0.03. 
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Method 


G/G(1CT 13 yr- 1 ) 


Reference 


Lunar Laser Ranging 


4±9 


295 


Binary Pulsar 1913 + 16 


40 ±50 


EH 


Helioseismology 


0±16 




Big Bang nucleosynthesis 


0±4 


EH 23 



Table 5: Constancy of the gravitational constant. For binary pulsar data, the bounds are dependent 
upon the theory of gravity in the strong-field regime and on neutron star equation of state. Big-bang 
nucleosynthesis bounds assume specific form for time dependence of G. 



Several observational constraints can be placed on G/G, one kind coming from bounding the 
present rate of variation, another from bounding a difference between the present value and a past 
value. The first type of bound typically comes from lunar laser-ranging measurements, planetary 
radar-ranging measurements, and pulsar timing data. The second type comes from studies of the 
evolution of the Sun, stars and the Earth, big-bang nucleosynthesis, and analyses of ancient eclipse 
data. Recent results are shown in Table El 

The best limits on a current G/G come from lunar laser ranging measurements (for earlier 
results see |9"51 12941 1193| 1. These have largely supplanted earlier bounds from ranging to the 
1976 Viking landers (see TEGP, 14.3 (c) |281| L which were limited by uncertain knowledge of the 
masses and orbits of asteroids. However, improvements in knowledge of the asteroid belt, combined 
with continuing radar observations of planets and spacecraft, notably the Mars Global Surveyor 
(1998-2003) and Mars Odyssey (2002 - present), may enable a bound on G/G at the level of a part 
in 10 13 per year. For an initial analysis along these lines, see |213) . It has been suggested that 
radar observations of the planned 2012 Bepi-Colombo Mercury orbiter mission over a two-year 
integration with 6 cm rms accuracy in range could yield A (G/G) < 10~ 13 yr -1 ; an eight-year 
mission could improve this by a factor 15 [1881 119| . 

Although bounds on G/G from solar-system measurements can be correctly obtained in a 
phenomenological manner through the simple expedient of replacing G by Go + Go(t — to) in 
Newton's equations of motion, the same does not hold true for pulsar and binary pulsar timing 
measurements. The reason is that, in theories of gravity that violate SEP, such as scalar-tensor 
theories, the "mass" and moment of inertia of a gravitationally bound body may vary with variation 
in G. Because neutron stars are highly relativistic, the fractional variation in these quantities 
can be comparable to AG/G, the precise variation depending both on the equation of state of 
neutron star matter and on the theory of gravity in the strong-field regime. The variation in the 
moment of inertia affects the spin rate of the pulsar, while the variation in the mass can affect the 
orbital period in a manner that can subtract from the direct effect of a variation in G, given by 
Pb/Pb = -2G/G HHJ. Thus, the bounds quoted in TableEJfor the binary pulsar PSR 1913+16 
and others |144| (see also are theory-dependent and must be treated as merely suggestive. 

In a similar manner, bounds from helioseismology and big-bang nucleosynthesis (BBN) assume 
a model for the evolution of G over the multi-billion year time spans involved. For example, the 
concordance of predictions for light elements produced around 3 minutes after the big bang with 
the abundances observed indicate that G then was within 20 percent of G today. Assuming a 
power-law variation of G ~ t~ a , then yields a bound on G/G today shown in Table El 
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3.7 Other tests of post-Newtonian gravity 



3.7.1 Search for gravitomagnetism 

According to GR, moving or rotating matter should produce a contribution to the gravitational field 
that is the analogue of the magnetic field of a moving charge or a magnetic dipole. In particular, 
one can view the goi part of the PPN metric (Box 2) as an analogue of the vector potential of 
electrodynamics. In a suitable gauge, and dropping the preferred-frame terms, it can be written 



At PN order, this contributes a Lorentz-type acceleration v x B s to the equation of motion, where 
the gravitomagnetic field B g is given by B g — V x (g oi e l ). 

Gravitomagnetism plays a role in a variety of measured relativistic effects involving moving 
material sources, such as the Earth-Moon system and binary pulsar systems. Nordtvedt 200 
1199] has argued that, if the gravitomagnetic potential H53[) were turned off, then there would be 
anomalous orbital effects in LLR and binary pulsar data. 

Rotation also produces a gravitomagnetic effect, since for a rotating body, V = — x J/r 3 , 
where J is the angular momentum of the body. The result is a "dragging of inertial frames" around 
the body, also called the Lense-Thirring effect. A consequence is a precession of a gyroscope's spin 
S according to 



where n is a unit radial vector, and r is the distance from the center of the body (TEGP 9.1 |281] \ 
The Relativity Gyroscope Experiment (Gravity Probe B or GPB) at Stanford University, in col- 
laboration with NASA and Lockheed-Martin Corporation 247 , recently completed a space mission 
to detect this frame-dragging or Lense-Thirring precession, along with the "geodetic" precession 
(see Sec. I3.7.2[l . A set of four superconducting- niobium-coated, spherical quartz gyroscopes were 
flown in a polar Earth orbit (642 km mean altitude, 0.0014 eccentricity), and the precessions of 
the gyroscopes relative to a distant guide star (HR 8703, IM Pegasi) were measured. For the given 
orbit, the predicted secular angular precession of the gyroscopes is in a direction perpendicular 
to the orbital plane at a rate ^(1 + 7 + jQ.%) x 41 x 10~ 3 arcsec/yr. The accuracy goal of the 
experiment is about 0.5 milliarcseconds per year. The spacecraft was launched on April 20, 2004, 
and the mission ended in September 2005, as scheduled, when the remaining liquid helium boiled 



It is too early to know whether the relativistic precessions were measured in the amount pre- 
dicted by GR, because an important calibration of the instrument exploits the effect of the aber- 
ration of starlight on the pointing of the on-board telescope toward the guide star, and completing 
this calibration required the full mission data set. In addition, part of the measured effect includes 
the motion of the guide star relative to distant inertial frames. This was measured before, during 
and after the mission separately by radio astronomers at Harvard/SAO and elsewhere using VLBI, 
and the results of those measurements were to be strictly embargoed until the GPB team has 
completed its analysis of the gyro data. Final results from the experiment are expected in 2006. 

Another way to look for frame-dragging is to measure the precession of orbital planes of bodies 
circling a rotating body. One implementation of this idea is to measure the relative precession, 
at about 31 milliarcseconds per year, of the line of nodes of a pair of laser-ranged geodynamics 
satellites (LAGEOS), ideally with supplementary inclination angles; the inclinations must be sup- 
plementary in order to cancel the dominant (126 degrees per year) nodal precession caused by the 
Earth's Newtonian gravitational multipole moments. Unfortunately, the two existing LAGEOS 
satellites are not in appropriately inclined orbits, and no concrete plans exist at present to launch 



90i = --(47 + 4 + a!)U 4 . 
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a third satellite in a supplementary orbit. Nevertheless, Ciufolini and Pavlis [Si] combined nodal 
precession data from LAGEOS I and II with improved models for the Earth's multipole moments 
provided by two recent orbiting geodesy satellites, Europe's CHAMP (Challenging Minisatellite 
Payload) and NASA's GRACE (Gravity Recovery and Climate Experiment), and reported a 5 - 
10 percent confirmation of GR. In earlier reports, Ciufolini et al. had reported tests at the the 20 
- 30 percent level, without the benefit of the GRACE/CHAMP data (SZIEH1EE1- Some authors 
stressed the importance of adequately assessing systematic errors in the LAGEOS data |227II134] . 

3.7.2 Geodetic precession 

A gyroscope moving through curved spacetime suffers a precession of its spin axis given by 



where v is the velocity of the gyroscope, and U is the Newtonian gravitational potential of the 
source (TEGP 9.1 |281p . The Earth-Moon system can be considered as a "gyroscope", with its 
axis perpendicular to the orbital plane. The predicted precession is about 2 arcseconds per century, 
an effect first calculated by de Sitter. This effect has been measured to about 0.6 percent using 
lunar laser ranging data [531 12941 l2*9"5] . 

For the GPB gyroscopes orbiting the Earth, the precession is 6.6 arcseconds per year. A goal 
of GPB is to measure this effect to 8 x 10~ 5 ; if achieved, this could bound the parameter 7 to a 
part in 10 4 , not competitive with the Cassini bound. 

3.7.3 Tests of post-Newtonian conservation laws 

Of the five "conservation law" PPN parameters £1, £2, £3, £4, and 0:3, only three, £2, C3 an d <%3, have 
been constrained directly with any precision; (1 is constrained indirectly through its appearance in 
the Nordtvedt effect parameter 77N, Eq. Ij51|l . There is strong theoretical evidence that (4, which is 
related to the gravity generated by fluid pressure, is not really an independent parameter - in any 
reasonable theory of gravity there should be a connection between the gravity produced by kinetic 
energy (pw 2 ), internal energy (pll), and pressure (p). From such considerations, there follows |275) 
the additional theoretical constraint 



A non-zero value for any of these parameters would result in a violation of conservation of 
momentum, or of Newton's third law in gravitating systems. An alternative statement of Newton's 
third law for gravitating systems is that the "active gravitational mass" , that is the mass that 
determines the gravitational potential exhibited by a body, should equal the "passive gravitational 
mass" , the mass that determines the force on a body in a gravitational field. Such an equality 
guarantees the equality of action and reaction and of conservation of momentum, at least in the 
Newtonian limit. 

A classic test of Newton's third law for gravitating systems was carried out in 1968 by Kreuzer, 
in which the gravitational attraction of fluorine and bromine were compared to a precision of 5 
parts in 10 5 . 

A remarkable planetary test was reported by Bartlett and van Buren [21]. They noted that 
current understanding of the structure of the Moon involves an iron-rich, aluminum-poor mantle 
whose center of mass is offset about 10 km from the center of mass of an aluminum-rich, iron- 
poor crust. The direction of offset is toward the Earth, about 14° to the east of the Earth-Moon 
line. Such a model accounts for the basaltic maria which face the Earth, and the aluminum-rich 
highlands on the Moon's far side, and for a 2 km offset between the observed center of mass 





(55) 
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and center of figure for the Moon. Because of this asymmetry, a violation of Newton's third 
law for aluminum and iron would result in a momentum non-conserving self-force on the Moon, 
whose component along the orbital direction would contribute to the secular acceleration of the 
lunar orbit. Improved knowledge of the lunar orbit through lunar laser ranging, and a better 
understanding of tidal effects in the Earth-Moon system (which also contribute to the secular 
acceleration) through satellite data, severely limit any anomalous secular acceleration, with the 
resulting limit 

(maK)ai - (m A /m P ) Fc 
(m A /m P ) Fc 

According to the PPN formalism, in a theory of gravity that violates conservation of momentum, 
but that obeys the constraint of Eq. (|56l) , the electrostatic binding energy E c of an atomic nucleus 
could make a contribution to the ratio of active to passive mass of the form 

m A = m P + 7j^(3E c . (58) 

The resulting limit on ( 3 from the lunar experiment is ( 3 < 1 x 1CT 8 (TEGP 9.2, 14.3 (d) [55T] ). 
Nordtvedt |204| has examined whether this bound could be improved by considering the asymmetric 
distribution of ocean water on Earth. 

Another consequence of a violation of conservation of momentum is a self-acceleration of the 
center of mass of a binary stellar system, given by 

1 mfiSm e 

a C M = -r(C2 +a 3 ) — T 5w2 np ' 59 

2 er a m (1 — e^yi z 

where 8m = m\ — rri2, a is the semi- major axis, and np is a unit vector directed from the center of 
mass to the point of periastron of mi (TEGP 9.3 281 ). A consequence of this acceleration would be 
non- vanishing values for cPP/dt 2 , where P denotes the period of any intrinsic process in the system 
(orbit, spectra, pulsar periods). The observed upper limit on d 2 P p /dt 2 of the binary pulsar PSR 
1913+16 places a strong constraint on such an effect, resulting in the bound \a 3 + (2] < 4 x 10~ 5 . 
Since 013 has already been constrained to be much less than this (Table 0}, we obtain a strong 
solitary bound on C2 < 4 x 10 -5 |2~5U] . 



< 4 x 10 
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3.8 Prospects for improved PPN parameter values 

A number of advanced experiments or space missions are under development or have been proposed 
which could lead to significant improvements in values of the PPN parameters, of 3 2 of the Sun, 
and of G/G. 

Lunar laser ranging at the Apache Point Observatory (APOLLO project) could improve bounds 
on the Nordvedt parameter to the level 3 x 10~ 5 and on G/G to better than 10 _13 yr _1 296 . 

The proposed 2012 ESA Bepi-Columbo Mercury orbiter, in a two-year experiment, with 6 cm 
range capability, could yield improvements in 7 to 3 x 10~ 5 , in (3 to 3 x 10 -4 , in ai to 10 , in 
G/G to 10 -13 yr _1 , and in J2 to 3 x 10 -8 . An eight-year mission could yield further improvements 
by factors of 2 - 5 in (3, ct\ and J2, and a further factor 15 in G/G |188llT§] . 

GAIA is a high- precision astrometric orbiting telescope (a successor to Hipparcos) , which could 
measure light-deflection and 7 to the 10 -6 level It is planned for launch by ESA in the 2011 
time frame. 

LATOR, Laser Astrometric Test of Relativity, is a concept for a NASA mission in which two 
microsatellites orbit the Sun on Earth-like orbits near superior conjunction, so that their lines of 
sight are close to the Sun. Using optical tracking and an optical interferometer on the International 
Space Station, it may be possible to measure the deflection of light with sufficient accuracy to bound 
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7 to a part in 10 and Ji to a part in 10 , and to measure the solar frame-dragging effect to one 
percent |5S0II35I] . 

Nordtvedt |203| has argued that "grand fits" of large solar system ranging data sets, including 
radar ranging to Mercury, Mars and satellites, and laser ranging to the Moon, could yield sub- 
stantially improved measurements of PPN parameters. A recent contribution in that direction 
is [5T3] . 

4 Strong Gravity and Gravitational Waves: A New Testing 
Ground 

4.1 Strong-field systems in general relativity 
4.1.1 Defining weak and strong gravity 

In the solar system, gravity is weak, in the sense that the Newtonian gravitational potential and 
related variables (E/(x, t) ~ v 2 ~ p/p ~ e) are much smaller than unity everywhere. This is the 
basis for the post-Newtonian expansion and for the "parametrized post-Newtonian" framework 
described in Sec. 13.21 "Strong-field" systems are those for which the simple 1PN approximation of 
the PPN framework is no longer appropriate. This can occur in a number of situations: 

• The system may contain strongly relativistic objects, such as neutron stars or black holes, near 
and inside which e ~ 1, and the post-Newtonian approximation breaks down. Nevertheless, 
under some circumstances, the orbital motion may be such that the interbody potential and 
orbital velocities still satisfy e < 1 so that a kind of post-Newtonian approximation for the 
orbital motion might work; however, the strong-field internal gravity of the bodies could 
(especially in alternative theories of gravity) leave imprints on the orbital motion. 

• The evolution of the system may be affected by the emission of gravitational radiation. The 
1PN approximation does not contain the effects of gravitational radiation back-reaction. In 
the expression for the metric given in Box [21 radiation back-reaction effects do not occur 
until 0(e 7 / 2 ) in goo, 0(e 3 ) in goi, and 0(e 5 / 2 ) in g^. Consequently, in order to describe such 
systems, one must carry out a solution of the equations substantially beyond 1PN order, 
sufficient to incorporate the leading radiation damping terms at 2.5PN order. In addition, 
the PPN metric described in Sec. 13.21 is valid in the near zone of the system, i.e. within 
one gravitational wavelength of the system's center of mass. As such it cannot describe the 
gravitational waves seen by a detector. 

• The system may be highly relativistic in its orbital motion, so that U ~ v 2 ~ 1 even for the 
interbody field and orbital velocity. Systems like this include the late stage of the inspiral 
of binary systems of neutron stars or black holes, driven by gravitational radiation damping, 
prior to a merger and collapse to a final stationary state. Binary inspiral is one of the leading 
candidate sources for detection by a world-wide network of laser interferometric gravitational 
wave observatories nearing completion. A proper description of such systems requires not 
only equations for the motion of the binary carried to extraordinarily high PN orders (at 
least 3.5PN), but also requires equations for the far-zone gravitational waveform measured 
at the detector, that are equally accurate to high PN orders beyond the leading "quadrupolc" 
approximation. 

Of course, some systems cannot be properly described by any post-Newtonian approximation 
because their behavior is fundamentally controlled by strong gravity. These include the imploding 
cores of supernovae, the final merger of two compact objects, the quasinormal-mode vibrations of 
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neutron stars and black holes, the structure of rapidly rotating neutron stars, and so on. Phe- 
nomena such as these must be analysed using different techniques. Chief among these is the 
full solution of Einstein's equations via numerical methods. This field of "numerical relativity" 
is a rapidly growing and maturing branch of gravitational physics, whose description is beyond 
the scope of this review (see [1661 126) for reviews) . Another is black hole perturbation theory 
(see [TBI da EU for reviews) . 

4.1.2 Compact bodies and the strong equivalence principle 

When dealing with the motion and gravitational-wave generation by orbiting bodies, one finds 
a remarkable simplification within GR. As long as the bodies are sufficiently well-separated that 
one can ignore tidal interactions and other effects that depend upon the finite extent of the bodies 
(such as their quadrupole and higher multipole moments) , then all aspects of their orbital behavior 
and gravitational wave generation can be characterized by just two parameters: mass and angular 
momentum. Whether their internal structure is highly relativistic, as in black holes or neutron 
stars, or non-relativistic as in the Earth and Sun, only the mass and angular momentum are needed. 
Furthermore, both quantities are measurable in principle by examining the external gravitational 
field of the bodies, and make no reference whatsoever to their interiors. 

Damour calls this the "effacement" of the bodies' internal structure. It is a consequence of 
the strong equivalence principle (SEP), described in Section 1.21 

General relativity satisfies SEP because it contains one and only one gravitational field, the 
spacetime metric g^ v . Consider the motion of a body in a binary system, whose size is small 
compared to the binary separation. Surround the body by a region that is large compared to 
the size of the body, yet small compared to the separation. Because of the general covariance of 
the theory, one can choose a freely-falling coordinate system which comoves with the body, whose 
spacetime metric takes the Minkowski form at its outer boundary (ignoring tidal effects generated 
by the companion). There is thus no evidence of the presence of the companion body, and the 
structure of the chosen body can be obtained using the field equations of GR in this coordinate 
system. Far from the chosen body, the metric is characterized by the mass and angular momentum 
(assuming that one ignores quadrupole and higher multipole moments of the body) as measured 
far from the body using orbiting test particles and gyroscopes. These asymptotically measured 
quantities are oblivious to the body's internal structure. A black hole of mass m and a planet of 
mass m would produce identical spacetimes in this outer region. 

The geometry of this region surrounding the one body must be matched to the geometry 
provided by the companion body. Einstein's equations provide consistency conditions for this 
matching that yield constraints on the motion of the bodies. These are the equations of motion. 
As a result the motion of two planets of mass and angular momentum mi, m2, Ji and J2 is 
identical to that of two black holes of the same mass and angular momentum (again, ignoring tidal 
effects). 

This effacement does not occur in an alternative gravitional theory like scalar-tensor gravity. 
There, in addition to the spacetime metric, a scalar field 4> is generated by the masses of the 
bodies, and controls the local value of the gravitational coupling constant (i.e. G is a function of 
<f) . Now, in the local frame surrounding one of the bodies in our binary system, while the metric 
can still be made Minkowskian far away, the scalar field will take on a value 4>q determined by the 
companion body. This can affect the value of G inside the chosen body, alter its internal structure 
(specifically its gravitational binding energy) and hence alter its mass. Effectively, each body can 
be characterized by several mass functions m^(i^), which depend on the value of the scalar field at 
its location, and several distinct masses come into play, such as inertial mass, gravitational mass, 
"radiation" mass, etc. The precise nature of the functions will depend on the body, specifically 
on its gravitational binding energy, and as a result, the motion and gravitational radiation may 
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depend on the internal structure of each body. For compact bodies such as neutron stars and black 
holes these internal structure effects could be large; for example, the gravitational binding energy 
of a neutron star can be 10 - 20 percent of its total mass. At 1PN order, the leading manifestation 
of this phenomenon is the Nordtvedt effect. 

This is how the study of orbiting systems containing compact objects provides strong-field 
tests of GR. Even though the strong-field nature of the bodies is effaced in GR, it is not in other 
theories, thus any result in agreement with the predictions of GR constitutes a kind of "null" test 
of strong-field gravity. 

4.2 Motion and gravitational radiation in general relativity 

The motion of bodies and the generation of gravitational radiation are long-standing problems 
that date back to the first years following the publication of GR, when Einstein calculated the 
gravitational radiation emitted by a laboratory-scale object using the linearized version of GR, 
and de Sitter calculated iV-body equations of motion for bodies in the 1PN approximation to GR. 
It has at times been controversial, with disputes over such issues as whether Einstein's equations 
alone imply equations of motion for bodies (Einstein, Infold and Hoffman demonstrated explicitly 
that they do, using a matching procedure similar to the one described above), whether gravitational 
waves are real or are artifacts of general covariance (Einstein waffled; Bondi and colleagues proved 
their reality rigorously in the 1950s), and even over algebraic errors (Einstein erred by a factor 
of 2 in his first radiation calculation; Eddington found the mistake). Shortly after the discovery 
of the binary pulsar PSR 1913+16 in 1974, questions were raised about the foundations of the 
"quadrupole formula" for gravitational radiation damping (and in some quarters, even about its 
quantitative validity). These questions were answered in part by theoretical work designed to 
shore up the foundations of the quadrupole approximation, and in part (perhaps mostly) by the 
agreement between the predictions of the quadrupole formula and the observed rate of damping of 
the pulsar's orbit (see Section 153)1 . Damour [72| gives a thorough historical and technical review 
of this subject up to 1986. 

The problem of motion and radiation in GR has received renewed interest since 1990, with 
proposals for construction of large-scale laser interferometric gravitational- wave observatories, such 
as the LIGO project in the US, VIRGO and GEO600 in Europe, and TAMA300 in Japan, and 
the realization that a leading candidate source of detectable waves would be the inspiral, driven 
by gravitational radiation damping, of a binary system of compact objects (neutron stars or black 
holes) 3, 257 . The analysis of signals from such systems will require theoretical predictions from 
GR that are extremely accurate, well beyond the leading-order prediction of Newtonian or even 
post-Newtonian gravity for the orbits, and well beyond the leading-order formulae for gravitational 
waves. 

This presented a major theoretical challenge: to calculate the motion and radiation of systems 
of compact objects to very high PN order, a formidable algebraic task, while addressing a number 
of issues of principle that have historically plagued this subject, sufficiently well to ensure that 
the results were physically meaningful. This challenge has been largely met, so that we may soon 
see a remarkable convergence between observational data and accurate predictions of gravitational 
theory that could provide new, strong-field tests of GR. 

Here we give a brief overview of the problem of motion and gravitational radiation in GR. 

4.3 Einstein's equations in "relaxed" form 

The Einstein equations G^ v — SirT^ are elegant and deceptively simple, showing geometry (in 
the form of the Einstein tensor G^u, which is a function of spacetime curvature) being generated 
by matter (in the form of the material stress-energy tensor T^). However, this is not the most 
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useful form for actual calculations. For post-Newtonian calculations, a far more useful form is the 
so-called "relaxed" Einstein equations: 

ah af3 = -16vrr Q/3 , (60) 

where □ = —d 2 /dt 2 + V 2 is the flat-spacetime wave operator, ft"' 3 is a "gravitational tensor 
potential" related to the deviation of the spacetime metric from its Minkowski form by the formula 
ff>P = 7y«/3 _ g) 1 / 2 g a P, g is the determinant of g a p, and a particular coordinate system has been 
specified by the de Dondcr or harmonic gauge condition dh a ^/dx 13 = (summation on repeated 
indices is assumed). This form of Einstein's equations bears a striking similarity to Maxwell's 
equations for the vector potential A a in Lorentz gauge: OA a = ~4irJ a , dA a /dx a = 0. There is a 
key difference, however: The source on the right hand side of Eq. I|60JI is given by the "effective" 
stress-energy pseudotensor 

T af3 = {-g)T a(} + (1677)"^^, (61) 

where A a @ is the non-linear "field" contribution given by terms quadratic (and higher) in h a @ and 
its derivatives (see |19(J) . Eqs. (20.20, 20.21) for formulae). In GR, the gravitational field itself 
generates gravity, a reflection of the nonlinearity of Einstein's equations, and in contrast to the 
linearity of Maxwell's equations. 

Equation 160f) is exact, and depends only on the assumption that spacetime can be covered by 
harmonic coordinates. It is called "relaxed" because it can be solved formally as a functional of 
source variables without specifying the motion of the source, in the form 

^(M)=4/ ^"-'*-^*'W , (02, 
Jc l x - x I 

where the integration is over the past flat-spacetime null cone C of the field point (t, x). The motion 
of the source is then determined either by the equation dr a ^/dx 13 = (which follows from the 
harmonic gauge condition), or from the usual covariant equation of motion T a @ .p = 0, where the 
subscript ; (3 denotes a covariant divergence. This formal solution can then be iterated in a slow 
motion (v < 1) weak-field (||ft Q ' 3 || < 1) approximation. One begins by substituting = into 
the source T a/3 in Eq. I|62|l . and solving for the first iterate /i" , and then repeating the procedure 
sufficiently many times to achieve a solution of the desired accuracy. For example, to obtain the 
1PN equations of motion, two iterations are needed (i.e. must be calculated); likewise, to 
obtain the leading gravitational waveform for a binary system, two iterations are needed. 

At the same time, just as in electromagnetism, the formal integral (|62() must be handled dif- 
ferently, depending on whether the field point is in the far zone or the near zone. For field points 
in the far zone or radiation zone, |x| > A > |x'| (A is the gravitational wavelength/27r), the field 
can be expanded in inverse powers of R = |x| in a multipole expansion, evaluated at the "retarded 
time" t — R. The leading term in 1/R is the gravitational waveform. For field points in the near 
zone or induction zone, |x| ~ |x'| < A, the field is expanded in powers of |x — x'| about the local 
time i, yielding instantaneous potentials that go into the equations of motion. 

However, because the source T af3 contains h a ^ itself, it is not confined to a compact region, 
but extends over all spacetime. As a result, there is a danger that the integrals involved in 
the various expansions will diverge or be ill-defined. This consequence of the non-linearity of 
Einstein's equations has bedeviled the subject of gravitational radiation for decades. Numerous 
approaches have been developed to try to handle this difficulty. The "post-Minkowskian" method 
of Blanchet, Damour and Iyer 1381 liTTA 1791 l4T)l 135] solves Einstein's equations by two different 
techniques, one in the near zone and one in the far zone, and uses the method of singular asymptotic 
matching to join the solutions in an overlap region. The method provides a natural "regularization" 
technique to control potentially divergent integrals (see [3Jy for a thorough review). The "Direct 
Integration of the Relaxed Einstein Equations" (DIRE) approach of Will, Wiseman and Pati [2911 
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209 retains Eq. l|62|l as the global solution, but splits the integration into one over the near zone 
and another over the far zone, and uses different integration variables to carry out the explicit 
integrals over the two zones. In the DIRE method, all integrals are finite and convergent. Itoh and 
Futamase have used an extension of the Einstein- Infeld-Hoffman matching approach combined with 
a specific method for taking a point-particle limit [135) . while Damour, Jaranowski and Schafer 
have pioneered an ADM Hamiltonian approach that focuses on the equations of motion |140l 11411 

EnnHaE] 

These methods assume from the outset that gravity is sufficiently weak that ||/i Q ^|| < 1 and 
harmonic coordinates exists everywhere, including inside the bodies. Thus, in order to apply the 
results to cases where the bodies may be neutron stars or black holes, one relies upon the strong 
equivalence principle to argue that, if tidal forces arc ignored, and equations are expressed in terms 
of masses and spins, one can simply extrapolate the results unchanged to the situation where the 
bodies are ultrarelativistic. While no general proof of this exists, it has been shown to be valid 
in specific circumstances, such as at 2PN order in the equations of motion, and for black holes 
moving in a Newtonian background field [72] . 

Methods such as these have resolved most of the issues that led to criticism of the foundations 
of gravitational radiation theory during the 1970s. 

4.4 Equations of motion and gravitational waveform 

Among the results of these approaches are formulae for the equations of motion and gravitational 
waveform of binary systems of compact objects, carried out to high orders in a PN expansion. 
Here we shall only state the key formulae that will be needed for this review. For example, the 
relative two-body equation of motion has the form 

a = — = -j {-n + Aip N + A 2PN + A 2 . 5 pn + A 3PN + A 3 . 5PN + ...}, (63) 
at r z 

where m = mi + to 2 is the total mass, r — |xi — x 2 |, v = vi — v 2 , and n = (xi — x 2 )/?\ The 
notation A„pn indicates that the term is 0(e n ) relative to the Newtonian term — n. Explicit 
and unambiguous formulae for non-spinning bodies through 3.5PN order have been calculated by 
various authors (see f° r a review). Here we quote only the first PN corrections and the leading 
radiation-reaction terms at 2.5PN order: 

A 1PN = J (4 + 2rj) ™ - (1 + 3f])v' 2 + ^r 2 j n + (4 - 2rj)r v, (64) 

A 2 . 5 PN = ~^]~ { (Sv 2 + 17— J rft - [3v 2 + 9— J vj , (65) 

where r\ = mxm%/ {mi + m 2 ) 2 . These terms are sufficient to analyse the orbit and evolution of the 
binary pulsar fSec. IS~T)l . For example, the 1PN terms are responsible for the periastron advance 
of an eccentric orbit, given by uj = Qirf^m/ a{l — e 2 ), where a and e are the semi-major axis and 
eccentricity, respectively, of the orbit, and /b is the orbital frequency, given to the needed order by 
Kepler's third law 2irf h = (to/o 3 ) 1 / 2 . 

Another product is a formula for the gravitational field far from the system, written schemati- 
cally in the form 

h i = —— - |q ■? + Q J 5PN + Q/pn + Ql.5PN + QJpn + Qi'sPN + •••}> (66) 

where R is the distance from the source, and the variables are to be evaluated at retarded time 
t — R. The leading term is the so-called quadrupole formula 

ft«(t,x) = |f«(t-fl), (67) 
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where F 3 is the quadrupole moment of the source, and overdots denote time derivatives. For a 
binary system this leads to 

Q l] = 2r](v l v j - mfi l h 3 /r). (68) 

For binary systems, explicit formulae for the waveform through 2PN order have been derived 
(see 02] for a ready-to- use presentation of the waveform for circular orbits; see for a full 
review) . 

Given the gravitational waveform, one can compute the rate at which energy is carried off by the 
radiation (schematically J hhdfl, the gravitational analog of the Poynting flux). The lowest-order 
quadrupole formula leads to the gravitational wave energy flux 

8 n TO 4 n , , 

E = -if — {I2v 2 -llf 2 ). 69 
15 

This has been extended to 3.5PN order beyond the quadrupole formula (see jHE] for a review). 
Formulae for fluxes of angular and linear momentum can also be derived. The 2.5PN radiation- 
reaction terms in the equation of motion tffirif) result in a damping of the orbital energy that precisely 
balances the energy flux l|69[) determined from the waveform. Averaged over one orbit, this results 
in a rate of increase of the binary's orbital frequency given by 

h = l ^f h {2nMKf/*F{el 

He) = ^-^(l+g^ + ge 4 ), (70) 

where M. is the so-called "chirp" mass, given by M. — jj 3 / 5 m. Notice that by making precise 
measurements of the phase $(t) = 2tt J f(t')dt' of either the orbit or the gravitational waves (for 
which / = 2/b for the dominant component) as a function of the frequency, one in effect measures 
the "chirp" mass of the system. 

These formalisms have also been generalized to include the leading effects of spin-orbit and 
spin-spin coupling between the bodies |146l 11451 1289| . 

Another approach to gravitational radiation is applicable to the special limit in which one mass 
is much smaller than the other. This is the method of black-hole perturbation theory. One begins 
with an exact background spacetime of a black hole, either the non-rotating Schwarzschild or the 
rotating Kerr solution, and perturbs it according to g^ v — gjQ + h^ u . The particle moves on a 
geodesic of the background spacetime, and a suitably defined source stress-energy tensor for the 
particle acts as a source for the gravitational perturbation and wave field h^ u . This method provides 
numerical results that are exact in v, as well as analytical results expressed as series in powers 
of v, both for non-rotating and for rotating black holes. For non-rotating holes, the analytical 
expansions have been carried to 5.5PN order, or e 5 ' 5 beyond the quadrupole approximation. All 
results of black hole perturbation agree precisely with the mi — ► limit of the PN results, up to 
the highest PN order where they can be compared (for reviews see j!89l I236] h 

4.5 Gravitational wave detection 

A gravitational wave detector can be modelled as a body of mass M at a distance L from a fiducial 
laboratory point, connected to the point by a spring of resonant frequency ojo and quality factor 
Q. From the equation of geodesic deviation, the infinitesimal displacement £ of the mass along the 
line of separation from its equilibrium position satisfies the equation of motion 

£ + 2^i + ufe = I (f+(9, 0, t/>)h+(t) + F x {9, </>, rP)h x (t)) , (71) 
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where F+(6,<fr,tp) and F x (9,(f>,if>) are "beam-pattern" factors, that depend on the direction of 
the source (9,(f>) and on a polarization angle ip, and h+(i) and h x (t) are gravitational waveforms 
corresponding to the two polarizations of the gravitational wave (for a review, see |256p . In a 
source coordinate system in which the x — y plane is the plane of the sky and the z-direction points 
toward the detector, these two modes are given by 

h+(t) = \{h^{t) - h y T y T (t)), h x (t) = h* v T (t), (72) 

where hl£ T represent transverse-traceless (TT) projections of the calculated waveform of Eq. (|fc>fc> |> . 
given by 

ht} T = h kl [{5 ik - N' l N k )(S jl - N j N l ) - -(6 ij - N i N j )(S kl - N k N 1 )], (73) 

where N 3 is a unit vector pointing toward the detector. The beam pattern factors depend on the 
orientation and nature of the detector. For a wave approaching along the laboratory z-direction, 
and for a mass whose location on the x — y plane makes an angle (f> with the x axis, the beam 
pattern factors are given by F + — cos 2(f> and F x = sin 2cf). For a resonant cylinder oriented along 
the laboratory z axis, and for source direction (9,(f)), they are given by F + = sin 2 9 cos 2ip, F x = 
sin 2 9sm2if> (the angle ip measures the relative orientation of the laboratory and source a;-axes). 
For a laser interferometer with one arm along the laboratory x-axis, the other along the y-axis, 
and with £ defined as the differential displacement along the two arms, the beam pattern functions 
are F + = i(l + cos 2 9) cos 2(f) cos 2ip — cos 9 sin 2(f) sin 2ip and F x = |(1 + cos 2 9) cos 2(f) sin 2if) + 
cos 9 sin 2(f) cos 2ifj. 

The waveforms h + (f) and h x (t) depend on the nature and evolution of the source. For example, 
for a binary system in a circular orbit, with an inclination i relative to the plane of the sky, and 
the x-axis oriented along the major axis of the projected orbit, the quadrupole approximation of 
Eq. (JHEl gives 

h+(t) = -^(2^/ b ) 2 / 3 (l+cos 2 z)cos2$ b (i), (74) 

h x (t) = -^-(27rA1/b) 2/3 2cosicos2$ b (i), (75) 
it 

where 3>b(^) = 2tt f* fb(t')dt' is the orbital phase. 



5 Stellar System Tests of Gravitational Theory 
5.1 The binary pulsar and general relativity 

The 1974 discovery of the binary pulsar B1913+16 by Joseph Taylor and Russell Hulse during a 
routine search for new pulsars provided the first possibility of probing new aspects of gravitational 
theory: the effects of strong relativistic internal gravitational fields on orbital dynamics, and the 
effects of gravitational radiation reaction. For reviews of the discovery, see the published Nobel 
Prize lectures by Hulse and Taylor [13311253] . For reviews of the current status of pulsars, including 
binary and millisecond pulsars, see |175II245] . 

The system consists of a pulsar of nominal period 59 ms in a close binary orbit with an as 
yet unseen companion. The orbital period is about 7.75 hours, and the eccentricity is 0.617. 
From detailed analyses of the arrival times of pulses (which amounts to an integrated version 
of the Doppler-shift methods used in spectroscopic binary systems), extremely accurate orbital 
and physical parameters for the system have been obtained (Table . Because the orbit is so 
close (« IRq) and because there is no evidence of an eclipse of the pulsar signal or of mass 
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Symbol 




Parameter 


(units) 


Value 


(i) "Physical" parameters 






Right Ascension 


a 


19 h 15 m 27. s 99999(2) 


Declination 


5 


16°06'27."4034(4) 


Pulsar Period 


P p (ms) 


59.0299983444181(5) 


Derivative of Period 




8.62713(8) x KT 18 


(ii) "Keplerian" parameters 






Projected semimajor axis 


dp sin i (s) 


2.341774(1) 


Eccentricity 


e 


0.6171338(4) 


Orbital Period 


P b (day) 


0.322997462727(5) 


Longitude of periastron 


(°) 


226.57518(4) 


Julian date of periastron 


T (MJD) 


46443.99588317(3) 


(iii) "Post-Keplerian" parameters 




Mean rate of periastron advance 


(u) (° yr^ 1 ) 


4.226595(5) 


Redshift /time dilation 


7' (ms) 


4.2919(8) 


Orbital period derivative 


p b (10- 12 ) 


-2.4184(9) 



Table 6: Parameters of the binary pulsar B1913+16. The numbers in parentheses denote errors 
in the last digit. Data taken from ]2(K \272\ . Note that 7' is not the same as the PPN parameter 7 
[see Eq. (7Zj)/. 



transfer from the companion, it is generally agreed that the companion is compact. Evolutionary 
arguments suggest that it is most likely a dead pulsar, while B1913+16 is a "recycled" pulsar. 
Thus the orbital motion is very clean, free from tidal or other complicating effects. Furthermore, 
the data acquisition is "clean" in the sense that by exploiting the intrinsic stability of the pulsar 
clock combined with the ability to maintain and transfer atomic time accurately using GPS, the 
observers can keep track of pulse time-of- arrival with an accuracy of 13 /is, despite extended gaps 
between observing sessions (including a several-year gap in the middle 1990s for an upgrade of the 
Arecibo radio telescope). The pulsar has shown no evidence of "glitches" in its pulse period. 

Three factors make this system an arena where relativistic celestial mechanics must be used: 
the relatively large size of relativistic effects [w rbit ~ (m/r) 1 / 2 » 10 -3 ], a factor of 10 larger than 
the corresponding values for solar-system orbits; the short orbital period, allowing secular effects 
to build up rapidly; and the cleanliness of the system, allowing accurate determinations of small 
effects. Because the orbital separation is large compared to the neutron stars' compact size, tidal 
effects can be ignored. Just as Newtonian gravity is used as a tool for measuring astrophysi- 
cal parameters of ordinary binary systems, so GR is used as a tool for measuring astrophysical 
parameters in the binary pulsar. 

The observational parameters that are obtained from a least-squares solution of the arrival-time 
data fall into three groups: (i) non-orbital parameters, such as the pulsar period and its rate of 
change (defined at a given epoch), and the position of the pulsar on the sky; (ii) five "Keplerian" 
parameters, most closely related to those appropriate for standard Newtonian binary systems, 
such as the eccentricity e, the orbital period Pb, and the semi-major axis of the pulsar projected 
along the line of sight, a p sin z; and (iii) five "post-Keplerian" parameters. The five post-Kcplerian 
parameters are: (d>), the average rate of periastron advance; 7', the amplitude of delays in arrival 
of pulses caused by the varying effects of the gravitational redshift and time dilation as the pulsar 
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moves in its elliptical orbit at varying distances from the companion and with varying speeds; Pb, 
the rate of change of orbital period, caused predominantly by gravitational radiation damping; and 
r and s = sini, respectively the "range" and "shape" of the Shapiro time delay of the pulsar signal 
as it propagates through the curved spacetime region near the companion, where i is the angle of 
inclination of the orbit relative to the plane of the sky. An additional 14 relativistic parameters 
are measurable in principle |9()j . 

In GR, the five post-Keplerian parameters can be related to the masses of the two bodies and 
to measured Keplerian parameters by the equations (TEGP 12.1, 14.6 (a) \2Sl] ) 

(cO) =6^/ b (2 7 rm/ b ) 2 / 3 (l-e 2 )- 1 , 
i = e(2 7 r/ b )- 1 (2^m/ b ) 2 / 3 (m 2 /m)(l + m 2 /m), 

P b = -(1927r/5)(27rX/ b ) 5 / 3 F(e), (76) 
s = sini, 
r = m 2 , 

where mi and m 2 denote the pulsar and companion masses, respectively. The formula for (uj) 
ignores possible non-relativistic contributions to the periastron shift, such as tidally or rotationally 
induced effects caused by the companion (for discussion of these effects, see TEGP 12.1 (c) [281J V 
The formula for P b includes only quadrupole gravitational radiation; it ignores other sources of 
energy loss, such as tidal dissipation (TEGP 12.1 (f) |281| V Notice that, by virtue of Kepler's 
third law, (27r/ b ) 2 = m/a 3 , (27rm/ b ) 2 / 3 = m/a ~ e, thus the first two post-Keplerian parameters 
can be seen as O(e), or 1PN corrections to the underlying variable, while the third is an 0(e 5 / 2 ), or 
2.5PN correction. The current observed values for the Keplerian and post-Keplerian parameters 
are shown in Table The parameters r and s are not separately measurable with interesting 
accuracy for B1913+16 because the orbit's 47° inclination does not lead to a substantial Shapiro 
delay. 

Because / b and e are separately measured parameters, the measurement of the three post- 
Keplerian parameters provide three constraints on the two unknown masses. The periastron shift 
measures the total mass of the system, P b measures the chirp mass, and 7' measures a complicated 
function of the masses. GR passes the test if it provides a consistent solution to these constraints, 
within the measurement errors. 

From the intersection of the (oj) and 7' constraints we obtain the values mi = 1.4414±O.OOO2Af 
and m 2 = 1.3867 ± 0.0002Af o . The third of Eqs. O then predicts the value P b = -2.40242 ± 
0.00002 x 10~ 12 . In order to compare the predicted value for P b with the observed value of Tabled 
it is necessary to take into account the small effect of a relative acceleration between the binary 
pulsar system and the solar system caused by the differential rotation of the galaxy. This effect 
was previously considered unimportant when P b was known only to 10 percent accuracy. Damour 
and Taylor |91| carried out a careful estimate of this effect using data on the location and proper 
motion of the pulsar, combined with the best information available on galactic rotation; the current 
value of this effect is Pg AL ~ -(0.0128 ± 0.0050) x 10" 12 . Subtracting this from the observed 
P b (Table EJ gives the corrected P^orr = _(2.4056 ± 0.0051) x 10~ 12 , which agrees with the 
prediction within the errors. In other words, 

pCORR 

b = 1.0013 ± 0.0021. (77) 

pGR 

The consistency among the measurements is displayed in Figure in which the regions allowed by 
the three most precise constraints have a single common overlap. Uncertainties in the parameters 
that go into the galactic correction are now the limiting factor in the accuracy of the test of 
gravitational damping. 
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A third way to display the agreement with GR is by comparing the observed phase of the orbit 
with a theoretical template phase as a function of time. If ft, varies slowly in time, then to first order 
in a Taylor expansion, the orbital phase is given by $t>(t) = 2nft,ot+TTft,ot 2 . The time of periastron 
passage tp is given by 3>(tp) = 2irN, where N is an integer, and consequently, the periastron 
time will not grow linearly with N. Thus the cumulative difference between periastron time tp 
and 7V//bo, the quantities actually measured in practice, should vary according to tp — N/ ft,o = 
—fwN 2 /2/^ — (/bo/2/bo)^ 2 - Figure d shows the results: the dots are the data points, while 
the curve is the predicted difference using the measured masses and the quadrupole formula for 

/bo eei. 

The consistency among the constraints provides a test of the assumption that the two bodies 
behave as "point" masses, without complicated tidal effects, obeying the general relativistic equa- 
tions of motion including gravitational radiation. It is also a test of strong gravity, in that the 
highly relativistic internal structure of the neutron stars does not influence their orbital motion, 
as predicted by the strong equivalence principle of GR. 




1.42 1.43 1.44 1.45 1.46 
MASS OF PULSAR (solar masses) 



Figure 6: Constraints on masses of the pulsar and its companion from data on B1913+16, assuming 
GR to be valid. The width of each strip in the plane reflects observational accuracy, shown as a 
percentage. An inset shows the three constraints on the full mass plane; the intersection region (a) 
has been magnified 400 times for the full figure. 

Recent observations |158l I271] indicate variations in the pulse profile, which suggests that the 
pulsar is undergoing geodetic precession on a 300-year timescale as it moves through the curved 
spacetime generated by its companion (Sec. I3.7.2I) . The amount is consistent with GR, assuming 
that the pulsar's spin is suitably misaligned with the orbital angular momentum. Unfortunately, 
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1975 1980 1985 1990 1995 2000 2005 
Year 

Figure 7: Plot of the cumulative shift of the periastron time from 1975-2005. The points are data, 
the curve is the GR prediction. The gap during the middle 1990s was caused by a closure of Arecibo 
for upgrading \272\ . 

the evidence suggests that the pulsar beam may precess out of our line of sight by 2025. 
5.2 A zoo of binary pulsars 

Nine relativistic binary neutron star systems with orbital periods less than a day are now known. 
While some are less interesting for testing relativity, some have yielded interesting tests, and 
others, notably the recently discovered "double pulsar" are likely to produce significant results in 
the future. Here we describe some of the more interesting or best studied cases; the parameters of 
the first four are listed in Tabled 

B 1534+12. This is a binary pulsar system in our galaxy |246ll245ll20| . Its pulses are significantly 
stronger and narrower than those of B1913+16, so timing measurements are more precise, 
reaching 3 /is accuracy. The orbital plane appears to be almost edge-on relative to the line 
of sight (i ~ 80°); as a result the Shapiro delay is substantial, and separate values of the 
parameters r and s have been obtained with interesting accuracy. Assuming GR, one infers 
that the two masses are mi = 1.335 ± 0.002M Q and m 2 = 1.344 ± 0.002M Q . The rate of 
orbit decay Pb agrees with GR to about 15 percent, but the precision is limited by the poorly 
known distance to the pulsar, which introduces a significant uncertainty into the subtraction 
of galactic acceleration. Independently of i\, measurement of the four other post-Keplerian 
parameters gives two tests of strong- field gravity in the non-radiative regime |254j . 

B2127+11C. This system appears to be a clone of the Hulse- Taylor binary pulsar, with very 
similar values for orbital period and eccentricity (see Tabled). The inferred total mass of 
the system is 2.706 ± 0. OHM©. But because the system is in the globular cluster M15 (NGC 
7078), it suffers Doppler shifts resulting from local accelerations, either by the mean cluster 
gravitational field or by nearby stars, that are more difficult to estimate than was the case 
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Parameter B1534+12 B2127+11C J1141-6545 J0737-3039(A,B) 

(i) "Keplerian" parameters 

a p sini (s) 3.7294626(8) 2.520(3) 1.85894(1) 1.41504(2)/1.513(3) 

e 0.2736767(1) 0.68141(2) 0.171876(2) 0.087779(5) 

P b (day) 0.420737299153(4) 0.335282052(6) 0.1976509587(3) 0.102251563(1) 

(ii) "Post-Keplerian" parameters 

(cJ) (° yr" 1 ) 1.755805(3) 4.457(12) 5.3084(9) 16.90(1) 

i (ms) 2.070(2) 4.67 0.72(3) 0.382(5) 

Pb (1(T 12 ) -0.137(3) -3.94 -0.43(10) -1.21(6) 

r(jus) 6.7(1.0) 6.2(5) 

s = sinz 0.975(7) 0.9995(4) 

Table 7: Parameters of other binary pulsars. References may be found in the text; for an online 
catalogue of pulsars with reasonably up-to-date parameters, see ]2(^ . 



with the galactic system B1913+16. This makes a separate, precision measurement of the 
relativistic contribution to Pb essentially impossible. 

J0737-3039A,B. This binary pulsar system, discovered in 2003 [SU], was already remarkable for 
its extraordinarily short orbital period (0.1 days) and large periastron advance (16.88° yr _1 ), 
but then the companion was also discovered to be a pulsar |176| . Because two projected 
semi-major axes can now be measured, one can obtain the mass ratio directly from the ratio 
of the two values of a p sin i, and thereby obtain the two masses by combining that ratio 
with the periastron advance, assuming GR. The results are toa = 1.337 ± 0.005 M© and 
tub = 1.250 ± 0.005 Mq, where A denotes the primary (first) pulsar. From these values, 
one finds that the orbit is nearly edge-on, with sin i = 0.9991, a value which is completely 
consistent with that inferred from the Shapiro delay parameter (Table 0|. In fact, the five 
measured post-Keplerian parameters plus the ratio of the projected semi-major axes give six 
constraints on the masses (assuming GR): all six overlap within their measurement errors. 
This system provides a unique opportunity for tight tests of strong-field and radiative effects 
in GR. Furthermore, it is likely that galactic proper motion effects will play a significantly 
smaller role in the interpretation of Pb measurements than they did in B1913+16. 

J1141-6545. This is a case where the companion is probably a white dwarf [521 EH- The mass 
of the pulsar and companion are 1.30 ± 0.02 and 0.986 ± 0.02 M©, respectively. P, has 
been measured to about 25 percent, consistent with the GR prediction. But because of the 
asymmetry in sensitivities (sns ~ 0.2, s^d ~ 10~ 4 ), there is the possibility, absent in the 
double neutron-star systems, to place a strong bound on scalar-tensor gravity (see Sec. 15. 4f) . 

J1756-2251. Discovered in 2004, this pulsar is in a binary system with a probable neutron star 
companion, with P b = 7.67 hr, e = 0.18 and lj = 2.585 ± 0.002 deg. yr" 1 [TUB] , 

J1906+0746. The discovery of this system was reported in late 2005 |174j . It is a young, 144-ms 
pulsar in a relativistic orbit with Pb = 3.98 hr, e = 0.085, and Co = 7.57 ± 0.03 deg. yr - . 

5.3 Binary pulsars and alternative theories 

Soon after the discovery of the binary pulsar it was widely hailed as a new testing ground for 
relativistic gravitational effects. As we have seen in the case of GR, in most respects, the system 



56 



has lived up to, indeed exceeded, the early expectations. 

In another respect, however, the system has only partially lived up to its promise, namely 
as a direct testing ground for alternative theories of gravity. The origin of this promise was the 
discovery that alternative theories of gravity generically predict the emission of dipole gravitational 
radiation from binary star systems. In GR, there is no dipole radiation because the "dipole 
moment" (center of mass) of isolated systems is uniform in time (conservation of momentum), 
and because the "inertial mass" that determines the dipole moment is the same as the mass 
that generates gravitational waves (SEP). In other theories, while the inertial dipole moment may 
remain uniform, the "gravity wave" dipole moment need not, because the mass that generates 
gravitational waves depends differently on the internal gravitational binding energy of each body 
than does the inertial mass (violation of SEP). Schematically, in a coordinate system in which the 
center of inertial mass is at the origin, so that mi.ixi + 771.1,2X2 = 0, the dipole part of the retarded 
gravitational field would be given by 



where v = vi — V2 and r\ and m are defined using inertial masses. In theories that violate 
SEP, the difference between gravitational-wave mass and inertial mass is a function of the internal 
gravitational binding energy of the bodies. This additional form of gravitational radiation damping 
could, at least in principle, be significantly stronger than the usual quadrupole damping, because 
it depends on fewer powers of the orbital velocity v, and it depends on the gravitational binding 
energy per unit mass of the bodies, which, for neutron stars, could be as large as 20 percent (see 
TEGP 10 |281) for further details). As one fulfillment of this promise, Will and Eardley worked 
out in detail the effects of dipole gravitational radiation in the bimetric theory of Rosen, and, 
when the first observation of the decrease of the orbital period was announced in 1979, the Rosen 
theory suffered a terminal blow. A wide class of alternative theories also fails the binary pulsar 
test because of dipole gravitational radiation (TEGP 12.3 |281p . 

On the other hand, the early observations of PSR 1913+16 already indicated that, in GR, the 
masses of the two bodies were nearly equal, so that, in theories of gravity that are in some sense 
"close" to GR, dipole gravitational radiation would not be a strong effect, because of the apparent 
symmetry of the system. The Rosen theory, and others like it, are not "close" to GR, except in their 
predictions for the weak-field, slow-motion regime of the solar system. When relativistic neutron 
stars are present, theories like these can predict strong effects on the motion of the bodies resulting 
from their internal highly relativistic gravitational structure (violations of SEP). As a consequence, 
the masses inferred from observations of the periastron shift and 7' may be significantly different 
from those inferred using GR, and may be different from each other, leading to strong dipole 
gravitational radiation damping. By contrast, the Brans-Dicke theory is "close" to GR, roughly 
speaking within 1/wbd of the predictions of the latter, for large values of the coupling constant 
wbd- Thus, despite the presence of dipole gravitational radiation, the binary pulsar provides at 
present only a weak test of Brans-Dicke theory, not competitive with solar-system tests. 

5.4 Binary pulsars and scalar-tensor gravity 

Making the usual assumption that both members of the system are neutron stars, and using 
the methods summarized in TEGP 10-12 |2"%T]. one can obtain formulas for the periastron shift, 
the gravitational redshift/second-order Doppler shift parameter, and the rate of change of orbital 
period, analogous to Eqs. I|76|l. These formulas depend on the masses of the two neutron stars, 
on their self-gravitational binding energy, represented by "sensitivities" s and k*, and on the 
Brans-Dicke coupling constant wbd- First, there is a modification of Kepler's third law, given by 




(78) 



2irf h = (Gm/a 3 ) 1 ^ 2 . 



(79) 
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Then, the predictions for (cl>), 7' and Pb are 

(u) = 67r/ b (2^m/b) 2/3 (l - e 2 )- 1 ^- 4 / 3 , (80) 

7' = e(2Trf h )- 1 {2irmf h ) 2 / 3 (m 2 /m)g- 1 / 3 (a* 2 +gm 2 /m + «^*), (81) 
P b = -(1927r/5)(27rX/b) 5/3 P'(e) - 4^(2^/ b )^ 2 G(e), (82) 
where M. = x 3 / 5 CJ~ 4 / 5 ?7 3 / 5 m, and, to first order in £ = (2 + cjbd) -1 , we have 



F'{e) = P(e) + + 3 r') 2 Qe 2 + ±e 4 ) (1 - e 2 )" 7 / 2 , 

G(e) = (1 - e 2 )- 5 /2 ^ + 1^ ; 



(83) 
(84) 



5 = ai-s 2 , Q = 1 - + s 2 - 2 Sl s 2 ) , V = G[l-^+^(s 1 + s 2 -2s 1 s 2 )}, (85) 



(86) 



c^ = l-£s 2 , r ?2 * = (l-2s 2 )£ ! x = ^ 

r = 1 - 2(mis 2 + m 2 si)/m , V = 1 — s\ — s 2 , (87) 
where P(e) is defined in Eq J7DJ). The quantities s a and k* are defined by 

/5(lnm a )\ „ /9(ln7 a 



■5'a 



V^lnG)^' a \d(lnG) /N 



and measure the "sensitivity" of the mass m a and moment of inertia I a of each body to changes 
in the scalar field (reflected in changes in G) for a fixed baryon number N (see TEGP 11, 12 and 
14.6 (c) |281| for further details). The quantity s a is related to the gravitational binding energy. 
These sensitivities will depend on the neutron-star equation of state. Notice how the violation 
of SEP in Brans-Dicke theory introduces complex structure-dependent effects in everything from 
the Newtonian limit (modification of the effective coupling constant in Kepler's third law) to 
gravitational radiation. In the limit £ — * 0, we recover GR, and all structure dependence disappears. 
The first term in P b (Eq. JB2J0 is the combined effect of quadrupole and monopole gravitational 
radiation, while the second term is the effect of dipole radiation. 

Unfortunately, because of the near equality of the neutron star masses in the binary pulsar, 
dipole radiation is suppressed, and the bounds obtained are not competitive with the Cassini bound 
on 7 |293j . except for those generalized scalar-tensor theories, with /3o < |77) . Bounds on the 
parameters ao and (3q from solar-system, binary-pulsar and gravitational wave observations (see 
Sees. I5.1l and lfi. 3(1 are found in [77) . 

Alternatively, a binary pulsar system with dissimilar objects, such as a white dwarf or black 
hole companion, would provide potentially more promising tests of dipole radiation. In this regard, 
the recently discovered binary pulsar J1141+6545, with an apparent white dwarf companion, may 
play an important role. Here one can treat %d ~ 10~ 4 as negligible. Then, from Eq. ((82(1 . it 
is straightforward to show that, if the timing reaches sufficient accuracy to determine Pb to an 
accuracy a in agreement with the prediction of GR, then the resulting lower bound on wbd would 
be 

w BD >4x 10 4 4 s /a. (89) 

Thus, for sns ~ 0.2, a 4 percent measurement would already compete with the Cassini bound (for 
further details, see |H5IITn3] V 
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6 Gravitational Wave Tests of Gravitational Theory 



6.1 Gravitational wave observatories 

Some time in the next decade, a new opportunity for testing relativistic gravity will be realized, 
when a worldwide network of kilometer-scale, laser interferometric gravitational wave observatories 
in the U.S. (LIGO project), Europe (VIRGO and GEO600 projects) and Japan (TAMA300 project) 
begins regular detection and analysis of gravitational-wave signals from astrophysical sources. 
These broad-band antennas will have the capability of detecting and measuring the gravitational 
waveforms from astronomical sources in a frequency band between about 10 Hz (the seismic noise 
cutoff) and 500 Hz (the photon counting noise cutoff), with a maximum sensitivity to strain at 
around 100 Hz of h ~ Al/l ~ 10~ 22 (rms), for the kilometer-scale LIGO/ VIRGO projects. The 
most promising source for detection and study of the gravitational wave signal is the "inspiralling 
compact binary" - a binary system of neutron stars or black holes (or one of each) in the hnal 
minutes of a death spiral leading to a violent merger. Such is the fate, for example, of the Hulse- 
Taylor binary pulsar B1913+16 in about 300 Myr, or the "double pulsar" J0737-3039 in about 85 
Myr. Given the expected sensitivity of the "advanced LIGO" (around 2010), which could see such 
sources out to many hundreds of megaparsecs, it has been estimated that from 40 to several hundred 
annual inspiral events could be detectable. Other sources, such as supernova core collapse events, 
instabilities in rapidly rotating newborn neutron stars, signals from non-axisymmetric pulsars, and 
a stochastic background of waves, may be detectable (for reviews, see |3l 1257) : for updates on the 
status of various projects, see |115llT7] '). 

A similar network of cryogenic resonant-mass gravitational antennas have been in operation 
for many years, albeit at lower levels of sensitivity (h ~ 10~ 19 ). While modest improvements 
in sensitivity may be expected in the future, these resonant detectors are not expected to be 
competitive with the large interferometers, unless new designs involving masses of spherical, or 
nearly spherical shape come to fruition. These systems are primarily sensitive to waves in relatively 
narrow bands about frequencies in the hundreds to thousands of Hz range |207l 11241 1341 1218) . 
although future improvements in sensitivity and increases in bandwidth may be possible |63| . 

In addition, plans are being developed for an orbiting laser interferometer space antenna (LISA 
for short). Such a system, consisting of three spacecraft orbiting the sun in a triangular formation 
separated from each other by five million kilometers, would be sensitive primarily in the very low 
frequency band between 10~ 4 and 10 _1 Hz, with peak strain sensitivity of order h ~ 10~ 23 |93) . 

In addition to opening a new astronomical window, the detailed observation of gravitational 
waves by such observatories may provide the means to test general relativistic predictions for 
the polarization and speed of the waves, for gravitational radiation damping and for strong-field 
gravity. 

6.2 Polarization of gravitational waves 

A laser interferometric or resonant bar gravitational wave detector whose scale is small compared 
to the gravitational wavelength measures the local components of a symmetric 3x3 tensor which 
is composed of the "electric" components of the Riemann curvature tensor, Roioj, via the equation 
of geodesic deviation, given, for a pair of freely falling particles by x 1 — —R^jX^, where x l denotes 
the spatial separation. In general there are six independent components, which can be expressed in 
terms of polarizations (modes with specific transformation properties under rotations and boosts). 
Three are transverse to the direction of propagation, with two representing quadrupolar deforma- 
tions and one representing a monopolar "breathing" deformation. Three modes are longitudinal, 
with one an axially symmetric stretching mode in the propagation direction, and one quadrupolar 
mode in each of the two orthogonal planes containing the propagation direction. Figure |H1 shows 
the displacements induced on a ring of freely falling test particles by each of these modes. General 
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relativity predicts only the first two transverse quadrupolar modes (a) and (b) independently of 
the source; these correspond to the waveforms h + and h x discussed earlier (note the cos 20 and 
sin 2(f> dependences of the displacements) . Massless scalar-tensor gravitational waves can in addi- 
tion contain the transverse breathing mode (c). In massive scalar-tensor theories, the longitudinal 
mode (d) can also be present, but is suppressed relative to (c) by a factor (A/A c ) 2 , where A is 
the wavelength of the radiation, and A c is the Compton wavelength of the massive scalar. More 
general metric theories predict additional longitudinal modes, up to the full complement of six 
(TEGP 10.2 fM) ). 

A suitable array of gravitational antennas could delineate or limit the number of modes present 
in a given wave. The strategy depends on whether or not the source direction is known. In general 
there are eight unknowns (six polarizations and two direction cosines), but only six measurables 
(RoiOj ) . If the direction can be established by either association of the waves with optical or other 
observations, or by time-of-flight measurements between separated detectors, then six suitably 
oriented detectors suffice to determine all six components. If the direction cannot be established, 
then the system is underdetermined, and no unique solution can be found. However, if one assumes 
that only transverse waves are present, then there are only three unknowns if the source direction 
is known, or five unknowns otherwise. Then the corresponding number (three or five) of detectors 
can determine the polarization. If distinct evidence were found of any mode other than the two 
transverse quadrupolar modes of GR, the result would be disastrous for GR. On the other hand, 
the absence of a breathing mode would not necessarily rule out scalar-tensor gravity, because the 
strength of that mode depends on the nature of the source. 

Some of the details of implementing such polarization observations have been worked out for 
arrays of resonant cylindrical, disk-shaped, spherical and truncated icosahcdral detectors (TEGP 
10.2 |281) . for recent reviews see [1701 [266 ); initial work has been done to assess whether the 
ground-based or space-based laser interferometers (or combinations of the two types) could perform 
interesting polarization measurements |267l 1491 11781 11181 1273) . Unfortunately for this purpose, the 
two LIGO observatories (in Washington and Louisiana states, respectively) have been constructed 
to have their respective arms as parallel as possible, apart from the curvature of the Earth; while 
this maximizes the joint sensitivity of the two detectors to gravitational waves, it minimizes their 
ability to detect two modes of polarization. 

6.3 Gravitational radiation back-reaction 

In the binary pulsar, a test of GR was made possible by measuring at least three relativistic 
effects that depended upon only two unknown masses. The evolution of the orbital phase under 
the damping effect of gravitational radiation played a crucial role. Another situation in which 
measurement of orbital phase can lead to tests of GR is that of the inspiralling compact binary 
system. The key differences are that here gravitational radiation itself is the detected signal, rather 
than radio pulses, and the phase evolution alone carries all the information. In the binary pulsar, 
the first derivative of the binary frequency /b, was measured; here the full nonlinear variation of 
/b as a function of time is measured. 

Broad-band laser interferometers are especially sensitive to the phase evolution of the gravi- 
tational waves, which carry the information about the orbital phase evolution. The analysis of 
gravitational wave data from such sources will involve some form of matched filtering of the noisy 
detector output against an ensemble of theoretical "template" waveforms which depend on the 
intrinsic parameters of the inspiralling binary, such as the component masses, spins, and so on, 
and on its inspiral evolution. How accurate must a template be in order to "match" the waveform 
from a given source (where by a match we mean maximizing the cross-correlation or the signal-to- 
noise ratio)? In the total accumulated phase of the wave detected in the sensitive bandwidth, the 
template must match the signal to a fraction of a cycle. For two inspiralling neutron stars, around 
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Figure 8: The six polarization modes for gravitational waves permitted in any metric theory of 
gravity. Shown is the displacement that each mode induces on a ring of test particles. The wave 
propagates in the +z direction. There is no displacement out of the plane of the picture. In (a), 
(b) and (c), the wave propagates out of the plane; in (d), (e), and (f), the wave propagates in 
the plane. In GR, only (a) and (b) are present; in massless scalar-tensor gravity, (c) may also be 
present. 
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16,000 cycles should be detected during the final few minutes of inspiral; this implies a phasing 
accuracy of 10~ 5 or better. Since v ~ 1/10 during the late inspiral, this means that correction 
terms in the phasing at the level of v 5 or higher are needed. More formal analyses confirm this 
intuition [fU EH EH ESI ■ 

Because it is a slow-motion system (v ~ 10 -3 ), the binary pulsar is sensitive only to the lowest- 
order effects of gravitational radiation as predicted by the quadrupole formula. Nevertheless, the 
first correction terms of order v and v 2 to the quadrupole formula were calculated as early as 



But for laser interferometric observations of gravitational waves, the bottom line is that, in order 
to measure the astrophysical parameters of the source and to test the properties of the gravitational 
waves, it is necessary to derive the gravitational waveform and the resulting radiation back-reaction 
on the orbit phasing at least to 2PN order beyond the quadrupole approximation, and preferably 
to 3PN order. 

For the special case of non-spinning bodies moving on quasi-circular orbits (i.e. circular apart 
from a slow inspiral), the evolution of the gravitational wave frequency / = 2/b through 2PN order 
has the form 



where r\ = 77*4 m 2 /m 2 . The first term is the quadrupole contribution (compare Eq. JZDJQ) the second 
term is the 1PN contribution, the third term, with the coefficient 47r, is the "tail" contribution, 
and the fourth term is the 2PN contribution, first reported jointly by Blanchet et al. |42l 1411 1291) . 
The 2.5PN, 3PN and 3.5PN contributions have also been calculated (see for a review). 

Similar expressions can be derived for the loss of angular momentum and linear momentum. 
Expressions for non-circular orbits have also been derived |122l 178) . These losses react back on the 
orbit to circularize it and cause it to inspiral. The result is that the orbital phase (and consequently 
the gravitational wave phase) evolves non-linearly with time. It is the sensitivity of the broad- 
band laser interferometric detectors to phase that makes the higher-order contributions to df jdt 
so observationally relevant. 

If the coefficients of each of the powers of / in Eq. l|9Ufl can be measured, then one again obtains 
more than two constraints on the two unknowns mj and m^, leading to the possibility to test GR. 
For example, Blanchet and Sathyaprakash [441 |4"5] have shown that, by observing a source with 
a sufficiently strong signal, an interesting test of the Air coefficient of the "tail" term could be 
performed. 

Another possibility involves gravitational waves from a small mass orbiting and inspiralling into 
a (possibly supermassive) spinning black hole. A general non-circular, non-equatorial orbit will 
precess around the hole, both in periastron and in orbital plane, leading to a complex gravitational 
waveform that carries information about the non-spherical, strong-field spacetime around the hole. 
According to GR, this spacetime must be the Kerr spacetime of a rotating black hole, uniquely 
specified by its mass and angular momentum, and consequently, observation of the waves could 
test this fundamental hypothesis of GR 232 214 . 

Thirdly, the dipole gravitational radiation predicted by scalar-tensor theories will result in 
a modification of the gravitational radiation back-reaction, and thereby of the phase evolution. 
Including only the leading quadrupole and dipole contributions, one obtains, in Brans-Dicke theory, 



1976 dU (see TEGP 10.3 |23T| 1. 




(90) 



f=^f(nMf)^ l + W)- 2/3 
5 



(91) 
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where M = (x 3 / 5 G~ 4 / 5 )r) 3 ^ m , and 6 is the coefficient of the dipole term, given by b = (5/48)(x _1 <? 4 / 3 )f<S 2 , 
where X) S are given by Eqs. JSJJ, and £ = 1/(2 + wbd)- Double neutron star systems are not 
promising because the small range of masses available near IAMq results in suppression of dipole 
radiation by symmetry For black holes, s = 0.5 identically, consequently double black-hole sys- 
tems turn out to be observationally identical in the two theories. Thus mixed systems involving a 
neutron star and a black hole are preferred. However, a number of analyses of the capabilities of 
both ground-based and space-based (LISA) observatories have shown that observing waves from 
neutron-star black-hole inspirals is not likely to bound scalar-tensor gravity at a level competitive 
with the Cassini bound or with future solar-system improvements 2831 11621 12371 l2*9~2l 1291 130) . 



6.4 Speed of gravitational waves 

According to GR, in the limit in which the wavelength of gravitational waves is small compared to 
the radius of curvature of the background spacetime, the waves propagate along null geodesies of the 
background spacetime, i.e. they have the same speed c as light (in this section, we do not set c = 1). 
In other theories, the speed could differ from c because of coupling of gravitation to "background" 
gravitational fields. For example, in the Rosen bimetric theory with a flat background metric 
77, gravitational waves follow null geodesies of r/, while light follows null geodesies of g (TEGP 

10.1 |2HU)- 

Another way in which the speed of gravitational waves could differ from c is if gravitation were 
propagated by a massive field (a massive graviton), in which case v g would be given by, in a local 
inertia! frame, 

2 2 4 

vi rate 

*-»--*-■ ^ w> 

where ra g and E are the graviton rest mass and energy, respectively. 

The simplest attempt to incorporate a massive graviton into general relativity in a ghost-free 
manner suffers from the so-called van Dam- Veltman-Zakharov (vDVZ) discontinuity [2631 [299 . 
Because of the 3 additional helicity states available to the massive spin-2 graviton, the limit of 
small graviton mass does not coincide with pure GR, and the predicted perihelion advance, for 
example, violates experiment. A model theory by Visser |265| attempts to circumvent the vDVZ 
problem by introducing a non-dynamical flat-background metric. This theory is truly continuous 
with GR in the limit of vanishing graviton mass; on the other hand, its observational implications 
have been only partially explored. Braneworld scenarios predict a tower or a continuum of massive 
gravitons, and may avoid the vDVZ discontinuity, although the full details are still a work in 
progress |MllfIH] . 

The most obvious way to test this is to compare the arrival times of a gravitational wave and 
an electromagnetic wave from the same event, e.g. a supernova. For a source at a distance D, the 
resulting value of the difference 1 — v g /c is 

V S = m-17 f 200 M P C \ f At 



'-f-»^'^JinJ' (93) 

where At = At a — (1 + Z)At e is the "time difference", where At a and At e are the differences in 
arrival time and emission time, respectively, of the two signals, and Z is the redshift of the source. 
In many cases, At e is unknown, so that the best one can do is employ an upper bound on At e 
based on observation or modelling. The result will then be a bound on 1 — v s /c. 

For a massive graviton, if the frequency of the gravitational waves is such that hf 3> ra g c 2 , 
where h is Planck's constant, then v g /c « 1 - ^(c/A g /) 2 , where A g = h/m g c is the graviton 
Compton wavelength, and the bound on 1 — v g /c can be converted to a bound on A g , given by 

A S > 3 x 10- km l^^ 1 -^) . (94) 



200 Mpc / J \fAt 
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The foregoing discussion assumes that the source emits both gravitational and electromagnetic 
radiation in detectable amounts, and that the relative time of emission can be established to 
sufficient accuracy, or can be shown to be sufficiently small. 

However, there is a situation in which a bound on the graviton mass can be set using gravi- 
tational radiation alone |285| . That is the case of the inspiralling compact binary. Because the 
frequency of the gravitational radiation sweeps from low frequency at the initial moment of obser- 
vation to higher frequency at the final moment, the speed of the gravitons emitted will vary, from 
lower speeds initially to higher speeds (closer to c) at the end. This will cause a distortion of the 
observed phasing of the waves and result in a shorter than expected overall time At a of passage of 
a given number of cycles. Furthermore, through the technique of matched filtering, the parameters 
of the compact binary can be measured accurately, (assuming that GR is a good approximation to 
the orbital evolution, even in the presence of a massive graviton), and thereby the emission time 
At e can be determined accurately. Roughly speaking, the "phase interval" fAt in Eq. (|94|l can 
be measured to an accuracy 1/p, where p is the signal-to-noise ratio. 

Thus one can estimate the bounds on A g achievable for various compact inspiral systems, and 
for various detectors. For stellar-mass inspiral (neutron stars or black holes) observed by the 
LIGO/VIRGO class of ground-based interferometers, D w 200 Mpc, / « fOO Hz, and fAt ~ 
ps f/fO. The result is A g > fO 13 km. For supermassive binary black holes (fO 4 to 10 7 M Q ) 
observed by the proposed laser interferometer space antenna (LISA), D « 3 Gpc, / ps 10~ 3 Hz, 
and fAt ~ p- 1 ps 1/1000. The result is A g > I0 17 km. 

A full noise analysis using proposed noise curves for the advanced LIGO and for LISA weakens 
these crude bounds by factors between two and 10 285 292, 23 ED]. For example, for the inspiral 
of two 10 6 Mq black holes with aligned spins at a distance of 3 Gpc observed by LISA, a bound 
of 2 x 10 16 km could be placed [2D]- Other possibilities include using binary pulsar data to 
bound modifications of gravitational radiation damping by a massive graviton |107) . and using 
LISA observations of the phasing of waves from compact white-dwarf binaries, eccentric galactic 
binaries, and eccentric inspiral binaries |71II143| . 

Bounds obtainable from gravitational radiation effects should be compared with the solid bound 
A g > 2.8 x 10 12 km, |251) derived from solar system dynamics, which limit the presence of a Yukawa 
modification of Newtonian gravity of the form 

V(r) = (GM/r) exp(-r/A g ), (95) 

and with the model-dependent bound A g > 6 x 10 19 km from consideration of galactic and cluster 
dynamics 265 . 

6.5 Strong-gravity tests 

One of the central difficulties of testing GR in the strong-field regime is the possibility of contami- 
nation by uncertain or complex physics. In the solar system, weak- field gravitational effects could 
in most cases be measured cleanly and separately from non-gravitational effects. The remarkable 
cleanliness of the binary pulsar permitted precise measurements of gravitational phenomena in a 
strong- field context. 

Unfortunately, nature is rarely so kind. Still, under suitable conditions, qualitative and even 
quantitative strong-field tests of GR could be carried out. 

One example is in cosmology. From a few seconds after the big bang until the present, the 
underlying physics of the universe is well understood, in terms of a Standard Model of a nearly 
spatially flat universe, 13.6 Byr old, dominated by dark matter and dark energy. Some alternative 
theories of gravity that are qualitatively different from GR fail to produce cosmologies that meet 
even the minimum requirements of agreeing qualitatively with big-bang nucleosynthesis (BBN) or 
the properties of the cosmic microwave background (TEGP 13.2 |281| 1. Others, such as Brans-Dicke 
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theory, are sufficiently close to GR (for large enough cjbd) that they conform to all cosmological 
observations, given the underlying uncertainties. The generalized scalar-tensor theories, however, 
could have small u> at early times, while evolving through the attractor mechanism to large u> today. 
One way to test such theories is through big-bang nucleosynthesis, since the abundances of the light 
elements produced when the temperature of the universe was about 1 MeV are sensitive to the rate 
of expansion at that epoch, which in turn depends on the strength of interaction between geometry 
and the scalar field. Because the universe is radiation-dominated at that epoch, uncertainties in 
the amount of cold dark matter or of the cosmological constant are unimportant. The nuclear 
reaction rates are reasonably well understood from laboratory experiments and theory, and the 
number of light neutrino families (3) conforms to evidence from particle accelerators. Thus, within 
modest uncertainties, one can assess the quantitative difference between the BBN predictions of 
GR and scalar-tensor gravity under strong-field conditions and compare with observations. For 
recent analyses, see |2351 1571 IrJUl 153] . 

Another example is the exploration of the spacetime near black holes and neutron stars via 
accreting matter. Studies of certain kinds of accretion known as advection-dominated accretion 
flow (ADAF) in low-luminosity binary X-ray sources may yield the signature of the black hole event 
horizon |186| . The spectrum of frequencies of quasi-periodic oscillations (QPO) from galactic black- 
hole binaries may permit measurement of the spins of the black holes |219) . Aspects of strong-field 
gravity and frame-dragging may be revealed in spectral shapes of iron fluorescence lines from the 
inner regions of accretion disks |226l 1225) . Because of uncertainties in the detailed models, the 
results to date of studies like these are suggestive at best, but the combination of higher- resolution 
observations and better modelling could lead to striking tests of strong-field predictions of GR. 

7 Conclusions 

We find that general relativity has held up under extensive experimental scrutiny. The question 
then arises, why bother to continue to test it? One reason is that gravity is a fundamental 
interaction of nature, and as such requires the most solid empirical underpinning we can provide. 
Another is that all attempts to quantize gravity and to unify it with the other forces suggest that 
the standard general relativity of Einstein is not likely to be the last word. Furthermore, the 
predictions of general relativity are fixed; the theory contains no adjustable constants so nothing 
can be changed. Thus every test of the theory is either a potentially deadly test or a possible probe 
for new physics. Although it is remarkable that this theory, born 90 years ago out of almost pure 
thought, has managed to survive every test, the possibility of finding a discrepancy will continue 
to drive experiments for years to come. 
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